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Abstract 

We propose an analytical approach to characterize undulator radiation near reso- 
nance, when the presence of the vacuum-pipe considerably affects radiation prop- 
erties. This is the case of the far-infrared undulator beamline at the Free-electron 
LASer (FEL) in Hamburg (FLASH), that will be capable of delivering pulses in 
the TeraHertz (THz) range. This undulator will allow pump-probe experiments 
where THz pulses are naturally synchronized to the VUV pulse from the FEL, as 
well as the development of novel electron-beam diagnostics techniques. Since the 
THz radiation diffraction-size exceeds the vacuum-chamber dimensions, charac- 
terization of infrared radiation must be performed accounting for the presence of 
a waveguide. We developed a theory of undulator radiation in a waveguide based 
on paraxial and resonance approximation. We solved the field equation with a 
tensor Green's function technique, and extracted figure of merits describing in a 
simple way the influence of the vacuum-pipe on the radiation pulse as a function 
of the problem parameters. Our theory, that makes consistent use of dimensionless 
analysis, allows treatment and physical understanding of many asymptotes of the 
parameter space, together with their region of applicability. 
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1 Introduction 



At the start of this century we have seen a revolution in synchrotron source 
intensities. This revolution stemmed from the technique of Free-Electron 
Lasers (FELs), combined with recent progress in accelerator technology, 
developed in connection with high-energy linear colliders. 

A new era of Synchrotron Radiation research has begun with first user 
experiments on VUV-FEL, based on Self- Amplified Spontaneous Emission 
(SASE) l^!,!?]. These results have been obtained at the TESLA (Tera Electron- 
volt Superconducting Linear Accelerator) Test Facility (TTF) at Deutsches 
Elektronen-SYnchrotron (DESY) at Hamburg, Germany. Radiation pulses 
with a wavelength of 98 nm, 40 fs duration and 1.5 GW peak power where 
used at TTF, phase 1 

This facility is now called FLASH (Free-Electron LAser in Hamburg), and 
operates as a user facility since August 2005 f3,^. Currently FLASH pro- 
duces GW-level, laser-like VUV (Vacuum Ultra- Violet) to EUV (Extreme 
Ultra- Violet) radiation pulses in the wavelength range between A = 13 nm 
and A = 50 nm, with a duration between 10 fs and 50 fs. Effective operation 
of an EEL in the VUV range requires high peak-current electron bunches (in 
the kA-range). Electrons are initially produced in a laser-driven RF-gun and 
subsequently compressed in magnetic chicanes to reach specifications. The 
accelerator complex at FLASH produces ultra-short bunches approaching 
sub-100 fs duration. 

It is natural to take advantage of these ultra-short bunches in order to 
provide coherent far-infrared (FIR) radiation. In fact, intense, coherent FIR 
radiation pulses can be produced from sub-100 fs electron bunches at a 
wavelength longer than, or comparable with the bunch length. For coher- 
ently radiating electrons, radiated energy is proportional to the square of 
the electron number, in contrast to the incoherent case when energy in the 
radiation pulse scales linearly with the number of electrons involved in the 
process. The result is an enhancement in the radiation intensity of many 
orders of magnitude. Reference llzj] describes a proposal for a FIR coherent 
source, integrated in the FLASH user facility. In [7] installation of an addi- 
tional electromagnetic undulator after the VUV EEL is proposed. The FIR 
source will use the spent electron beam coming from the EEL process, and 
will allow to significantly extend the scientific potential of FLASH without 
interfering with the main (VUV) option in the FLASH operation. 

Coherent FIR pulses will be intrinsically synchronized with the VUV pulses. 
A first, natural application of this kind of photon beams is for pump-probe 
experiments. Another application concerns electron beam diagnostics. The 
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femtosecond time-scale of electron bunches at FLASH is beyond the range 
of standard electronic-display instrumentation, and the development of 
non-destructive methods for the measurement of longitudinal beam cur- 
rent distribution is undoubtedly a challenging problem. In [8] a diagnostics 
technique is proposed, based on the measurement of FIR coherent radiation 
from electron bunches passing through an undulator. This technique is non- 
destructive, and characterization of bunches with strongly non-Gaussian 
shapes (as in the FLASH case) is possible 



FLASH will soon be coupled with a FIR electromagnetic undulator [11]. 
At the time of writing, the project is entering into the realization phase. 
Funding has been secured, and the undulator has been built and delivered 
|12|1 . Undulator radiation around the fundamental harmonic will always be 
used. The fundamental harmonic will be tuned by adjusting the magnetic 
field strength. The wavelength range (A = 60 -r 200 [im) provided by this 
powerful radiation source (up to 10 MW peak power) will overlap with 
a large part of the THz-gap, extending between 60 /.tm and 600 /^m. This 
will allow both applications for pump-probe experiments combining FIR 
and VUV radiation [13], and for non-destructive electron beam diagnostics. 
In fact, measurements of the modulus of the electron bunch form factor 
in the spectral range A = 10 200 /.im is sufficient to provide a precise 
reconstruction of the electron bunch profile |10|] . 



In the case of the FIR undulator beamline at FLASH the electron beam geo- 
metrical emittance is much smaller than the radiation diffraction size. This 
means that, as pertains the characterization of the THz pulses, the electron 
beam can safely be modelled as a filament beam. Computer codes like SRW 
1 14 ] and SPECTRA [jisl] can be used, in the space-frequency domain, to study 
undulator radiation (UR) from filament electron beams up to a wavelength 
when the influence of the vacuum chamber is negligible. Alternatively, an 
analytical formalism for describing near-zone Synchrotron Radiation (SR) 
fields from undulators in terms of Fourier optics, was developed in [16], 
also working in the free-space limit and in the space-frequency domain. 
These methods help designing beamlines and experiments. However, in the 
case of the FIR undulator beamline at FLASH, wavelengths in the order of 
200 [im and a 4 m-long undulator yield a radiation diffraction size of order of 
a centimeter. This rough estimate indicates that vacuum chamber effects are 
expected to play an important role, since a circular vacuum chamber with 
radius R = 1.8 cm is foreseen. In this case, conventional computer codes and 
analytical methods fail to predict the correct radiation characteristics. 



Summing up, in view of the practical application to the infrared undula- 
tor line at FLASH, there is a need to develop a comprehensive theory of 
undulator radiation in the presence of a waveguide. 
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Optimization of the radiation transport system calls for a precise character- 
ization of THz pulses along the photon beamline. In the present work we 
focus on the characterization of undulator radiation in presence of a waveg- 
uide, extending the study presented in [16] to include the influence of a 
vacuum chamber in our consideration. Our work is, therefore, propaedeu- 
tical to the problem of optimizing the radiation transport system after the 
undulator, but it does not directly deal with it. 

The task that one has to solve differs from the free-space case only in the 
formulation of boundary conditions. As we will discuss in Section 3, the 
paraxial approximation applies as in free-space but, on a perfectly conduc- 
tive boundary, the electric field must be orthogonal to the pipe surface. As in 
the free-space case one can use a Green's function approach to solve the field 
equations. The presence of different boundary conditions complicates the 
solution of the paraxial equation for the field, which can anyway be found 
explicitly, by accounting for the tensorial nature of the Green's function. 
Our consideration is quite general, and can be applied to UR sources as well 
as to other long-wavelength radiation setups, like edge-radiation setups. 
In both cases the paraxial approximation can be used. In the UR case, the 
resonance approximation can be exploited too, in addition to the paraxial 
approximation. Our theory is developed under both paraxial and resonance 
approximation. We thus consider a large number of undulator periods and 
a frequency range of interest close to the fundamental harmonic, where the 
free-space field exhibits horizontal polarization (for undulator field in the 
vertical direction) and azimuthal symmetry. The simultaneous application 
of paraxial and resonance approximation makes the UR case richer and 
more difficult to be studied from a theoretical viewpoint when compared 
with the edge-radiation case. In particular, the region of applicability of both 
approximations must be discussed and specified. 

In this paper we will focus on UR only, leaving the study of the edge- 
radiation case to a future, dedicated publication. In free-space and under 
resonance approximation UR is horizontally polarized. This is actually a 
replica of the undulator polarization properties. Moreover, the field exhibits 
azimuthal symmetry. These properties are lost when metallic boundaries 
are introduced. It is worth mentioning that application of similarity tech- 
niques helps to give a clear physical interpretation of numerical results. In 
this paper we continue to use these techniques, that we already applied to 
SR theory in free-space (see e.g. fl^). Despite the fact that equations for the 
undulator source are significantly complicated by the presence of a waveg- 
uide, we find that these complications result in the appearance of a single 
dimensionless extra-parameter, namely the waveguide diffraction parame- 
ter Q = R^/(ALi„), where A = A/{2n) is the reduced radiation wavelength. 
La, is the undulator length and, as already said, R is the waveguide radius. 
The physical interpretation of Q is the squared ratio between the waveguide 



5 



radius and the radiation diffraction size of UR in free-space. 



Since we are practically interested in the FIR undulator beamline at FLASH 
we put particular emphasis on planar undulators in the presence of a pipe 
with circularly symmetric cross-section. An explicit expression for the field 
is calculated as a superposition of Transverse Electric (TE) and Transverse 
Magnetic (TM) modes. Some figure of merit should be extracted from the full 
information carried by the expression for the field about how the metallic 
pipe influences radiation properties. We separately studied, for horizontal 
and vertical polarization components, two-dimensional intensity distribu- 
tions on a transverse plane at arbitrary distance from the undulator, for 
different choices of the problem parameter. Also, we analyzed the total 
power as a function of the waveguide diffraction parameter at perfect res- 
onance. Conversely, once the waveguide diffraction parameter is fixed, one 
can investigate how the total power changes as a function of the detuning 
from resonance. Finally, a comparison between the magnitude of the hori- 
zontally and vertically polarized fields is also proposed as a measure of the 
waveguide influence. 



To the best of our knowledge, only a few articles |17l, Il8l, |19|1 deal with the 
problem of radiation from a wiggled electron in a waveguide. This problem 
is also discussed at advanced textbook level in [20]. In reference [19] one 
may find the following words: "Motz and Nakamura [17] (...) considered 
an infinitely long wiggler. As a consequence, the outcome did not possess 
a realistic bandwidth. (...). In a (...) article by Haus and Islam [18] (...) it is 
shown that in the limit of a highly over-moded guide, a result similar to 
the free-space expression is produced". Reference [18] considers a rectan- 
gular waveguide, and focuses on "similarities between the emission into 
the free and the bounded space" (cited from [19]). However, we should say 
that reference |[l8|] solves the equations for the field in the rectangular case 
in all generality, and with no restrictions on the undulator parameter K. 
Reference [|19|] deals with differences with respect to free-space emission, in 
the case of a planar waveguide and planar undulator with small undulator 
parameter K «: 1. These restrictions limit the practical scope of that work. 
In particular, extending the theory in [19] from the case of a planar waveg- 
uide to other geometries is not straightforward, since the explicit expression 



for the paraxial Green's function given in |19l] is only valid for the planar- 
waveguide case. The planar-waveguide case has clear advantages from an 
educational viewpoint, as it allows to reduce complexities to a minimum. In 
particular, no restriction is made on the wiggling amplitude of the electron 
motion. In fact, the electron oscillates in the unbounded region between the 
parallel conducting plates. A transparent physical picture arises, in terms of 
reflection on the two metallic plates. However, such transparency comes at 
the cost of a limited r egio n of applicability of the theory in practical cases of 



interest. In reference ||19[1 , the resonance approximation is exploited, and a 
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spectral region for large detuning from the first harmonic can be considered 
due to simplifications intrinsic in the planar-waveguide geometry. In the 
case of the FLASH infrared undulator a circular waveguide is going to be 
used. The FLASH infrared undulator will operate for large values of K, and 
radiation will be used near the resonance with the first harmonic. We devel- 
oped a theory to deal with this situation by restricting our attention around 
resonance. Thus, we eliminated restrictions for the waveguide geometry 
(and for the undulator parameter K), but we introduced a condition about 
the spectral range of interest, coinciding with the spectral range of interest at 
the FLASH infrared undulator line. The wiggling amplitude of the electron 
in the undulator is taken to be small with respect to the dimension of the 
waveguide. This greatly simplifies analytical calculations, and describes our 
practical case of interest. As it will be seen in Section 5, within the region 
of parameter space where these conditions apply, one may consider both 
cases when the waveguide influence is weak (up to the free-space limit) 
or strong. Finally, the case of radiation of an electron in a helical wiggler 
with a circular waveguide is discussed in Section 9.3 of reference [20]. TM 
modes are neglected (see Appendix A for details). Excluding TM modes 
from consideration is held by us to be a misconception. In contrast to this, 
we will demonstrate throughout the text that the wiggler-induced motion 
couples with both TE and TM modes in a circular (cylindrical) waveguide. 

Our work is organized as follows. Besides this Introduction, in the next 
Section 2 we review basic theory of undulator radiation in free-space. Such 
a review is necessary, because delicate physical assumptions used in the 
free-space case continue to be exploited when a waveguide is present as 
well. In the following Section 3 we pose our problem, discussing Maxwell's 
equations for a single electron moving in the presence of metallic boundaries 
and within the paraxial approximation. We derive a closed expression for 
the field with the help of a tensor Green's function technique. At this stage 
the geometry of the vacuum pipe is generic, as well as the trajectory of 
the electron. In Section 4 we fix the vacuum pipe geometry, specializing our 
equations to the case of a vacuum pipe with circular cross-section. We verify 
the correctness of our results by deriving the free-space limit, studied in [16]. 
Then, in the following Sections 5 we consider the cases of a planar undulator, 
while we refer the interested reader to Appendix A for the case of a helical 
undulator. Results obtained in Section 5 for the planar undulator case are 
the main results in our paper, as they allow a complete characterization 
of the field in the FIR undulator beamline at FLASH. As said before, such 
characterization is indispensable for any further analysis of the radiation 
transport system. We analyze our results further in Section 6, where we 
propose a study of figure of merits, that will help designers and beamline 
scientists to estimate the influence of the vacuum pipe on the radiation 
characteristics. In Section 7 we discuss the influence of wall-resistance on 
our findings, demonstrating their relevance. Finally, in Section 8, we come 
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to conclusions. 



2 Undulator radiation in free-space 

2.2 Green's function technique in free-space 

For any SR setup we can represent the electric field in time domain j^(r, f) 
as a time-dependent function of an observation point located at position 
r = rj_+zez = xCx + yey + zcz, where e^, Cy and are defined as (dimen- 
sionless) unit vectors along horizontal, vertical and longitudinal direction 

in a given reference frame. For monochromatic waves of angular frequency 

_^ ^ 

CO, the wave amplitude has the form E{z,r^,t) = £(z,rL)exp[-/a;f] + C.C., 

■± 

where "C.C." indicates the complex conjugate of the preceding term and £ 

■± 

describes the variation of the wave amplitude. The vector £ actually rep- 
resents the amplitude of the electric field in the space-frequency domain. 

Accounting for electromagnetic sources, i.e. in a region of space where cur- 
rent and charge densities are present, the field in the space-frequency obeys 

-± ^ -> -± 

Helmholtz's equation, i.e. c^V^£ -I- aP-E - 4nc^Vp - Anicoj, where p{r, co) and 
-± 

j{r, co) are the Fourier transforms of the charge density, p{r, t), and of the 
current density, j{r, t), while c is the speed of light. We will consider a single 
electron. Using the Dirac delta distribution, we can write p{r, t) = —eb{r—r'{t)) 

and t) - v{t)p{r, t), where {-e) is the negative electron charge, r'{t) and v{t) 
are, respectively, the position and the velocity of the particle at a given time 
f, and Vz its longitudinal velocity. After calculation of the Fourier transform 
of these quantities and substitution into Helmholtz's equation we obtain 



lico 4ne 
H — I £ = — — exp 



ICO 



^i?(z)-v|6(r-l-r\(z)),(l) 



s(z) being the curvilinear abscissa measured along the electron trajectory, 
where we conventionally set s(0) = 0. When the longitudinal velocity of 
the electron, Vz, is close to the speed of light c (i.e. yl » 1, where yz(z) = 
(1 - vl/c^)~^^^ is the longitudinal Lorentz factor), the Fourier components of 
the source are almost synchronized with the electromagnetic wave travelling 
at the speed of light. In this case the phase co{s{z)/v-z/c) is a slow function of 
z compared to the wavelength. For example, in the particular case of motion 
on a straight section, one has s(z) = z/Vz, so that co{s{z)/v - z/c) - coz/{2ylc), 
and if ^ 1 such phase grows slowly in z with respect to the wavelength. 
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For a more generic motion, one similarly obtains: 



's{Z2) - s(Zi) 




(2) 



Mathematically, the phase in Eq. (2) enters in the Green's function solution 
of Eq. (1) as a factor in the integrand. As we integrate along z', the factor 
a){s{z')/v - z' jc) leads to an oscillatory behavior of the integrand over a 
certain integration range in z'. Such range can be identified with the value 
of Z2 - Zl for which the right hand side of Eq. (2) is of order unity, and 
it is naturally defined as the radiation formation length of the system 
at frequency co. Of course there exists some freedom in the choice of such 
definition: "order of unity" is not a precise number, and reflects the fact that 
there is no abrupt threshold between "oscillatory" and "non-oscillatory" 
behavior of the integrand in the solution of Eq. (1). In the following we 
define the formation length Ly^ as the interval Z2 - Zi such that the right hand 
side of Eq. (2) is strictly equal to unity. It is easy to see by inspection of 
Eq. (2) that if is sensibly smaller than c (but still of order c), i.e. ~ c 
but l/y^ ~ 1/ then L/ ~ A, where we introduced the reduced wavelength 
A = AIQji), and A = Incjco. On the contrary, when Vz is very close to c, i.e. 
llyl^l, the right hand side of Eq. (2) is of order unity for = Z2 - Zi » A. 
When the radiation formation length is much longer than A, the electric field 

envelope = Ej_ exp [-icoz/c] does not vary much along z on the scale of 
A, that is I dzEx,y co/c \ Ex,y |. Therefore, the second order derivative with 
respect to z in the operator on the left hand side of Eq. (1) is negligible 
with respect to the first order derivative with respect to z. As a result, Eq. 
(1) can be simplified as 



E±{z,fj_) 



= f{z,rj_) . 



(3) 



The differential operator D in Eq. (3) is defined by 



where V^^ is the Laplacian operator over transverse cartesian coordinates. 

The vector /(z, r_i) is specified by the trajectory of the source electron, r\{z), 
and is written as 
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r* 4:716 

f= — exp 



z 

I 



dz- 



co 



lyl{z)c 



6(rL -r\(z)) 



(5) 



Here we considered transverse components of £ only and we substituted 
Vz{z) with c, based on the fact that l/y^ ^ 1- Eq. (3) is Maxwell's equation 
in paraxial approximation. 

In the following, with some abuse of language, we will refer to the slowly 
varying envelope of Ex, i.e. to E^, simply as the field. 

The Green's function for Eq. (3), namely the solution corresponding to a 
unit point source, depends on boundary conditions. In free-space it must 
obey 



G(z - z'; rx - r'x) = 



47i(z - z') 



exp 



r', 



ICO- 



2c{z - z') 



(6) 



assuming z - z' > 0. When z - z' < the paraxial approximation does 
not hold, and the paraxial wave equation, Eq. (3), should be substituted, 
in the space-frequency domain, by the more general Helmholtz's equation. 
However, the radiation formation length for z - z' < is very short with 
respect to the case z - z' > 0, i.e. there is effectively no radiation for observer 
positions z-z' < 0. As a result, in this paper we will consider only z-z' > 0. 
It follows that the observer is located downstream of the sources. This leads 
to the final result: 



c- J z-z' 



v(z') rj_ - r'x(z') 
c z-z' 



xexp 



tco 



r± - r 



2c{z 



z') J lyl{z)c 



(7) 



Eq. (7) is valid at any observation position z such that the paraxial approx- 
imation is valid, i.e. up to distances between the observer and the elec- 
tromagnetic sources comparable with the radiation wavelength. One may 
recognize two terms in Eq. (7). The first in v(z') can be traced back to the 
current term in Eq. (1), while the second, in rx - ^'±(z')/ corresponds to the 
gradient term in Eq. (1). 

The far zone is defined as the region of observation where the direction from 
any trajectory point to the observer can be considered constant. In the far 
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zone, Eq. (7) reduces to 



xexp 



2' 

J 2yl{z)c 2c V ^ 

. 



(8) 



where 6 = r±/z defines the observation direction, and 6 = \d\. 

An ultra-relativistic electron radiates non-negligibly, in the far-zone, up 
to observation angles of order 0^ < A/Lf (see [16]). Moreover, for ultra- 
relativistic electrons, L^^ » A. The angular region 0^ < A/L^ ^ 1 formally 
coincides with the region of applicability of the paraxial approximation. 
However, since radiation is negligible elsewhere, it follows that paraxial 
approximation can be applied to describe radiation from an ultra-relativistic 
electron at any observation angle of interest. 

If we now decompose the electric field distribution in the far zone as a 

superposition of plane waves (angular spectrum), the angle of propagation 

of each plane wave is represented by the ratio between the transverse wave 

-» -> 
vector kj^ and the longitudinal wave number kz, that is kj_/kz. Note that in 
-> 

free-space, k_L and kz are allowed to vary continuously across the reciprocal 
space. Non-negligible plane wave components of the angular spectrum are 
those seen at angles 6^ < A/Lf « 1, i.e. propagating at angles k]_/kl « 1. 
Now, since kz - k = 1/A this may also be stated by requiring k]_c^/a)^ <c 1. 
This result will be useful in Section 3. 



2.2 Undulator radiation in paraxial approximation. Far zone 



Let us apply the method outlined in Section 2.1 to the case of planar undula- 
tor. In general, the term undulator radiation field means only a part of the total 
field seen by an observer from a realistic setup, because one should account 
for contributions from the entire trajectory of the particle. We may follow 
any textbook like [26] in deriving well-known relations in the far-zone. For 
the electron transverse velocity we assume 



c6s 

V:ciz') = -cds sm{ku,z') = --^ [exp{ik-,,z') - exp{-ik^z')] . (9) 
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Here kj^ = In I A^,, and A^; is the undulator period. Moreover, 0s = i^ly, where 
K is the deflection parameter defined as 



K = 



(10) 



Me being the electron mass at rest and Hy, being the maximal magnetic field 
of the undulator on axis. 

The longitudinal Lorentz factor yz{z) is a function of the position down the 
undulator, so that the phase in Eq. (2) calculated at Zi = and Z2 = z gives: 



z 

/ 



CO CO cod . 

-az = z — -z - — — sm(2K„,z) , 



2cyl{z) Icyl 8k^c 



(11) 



where the average longitudinal Lorentz factor yz is defined as 



fz = 



y 



(12) 



We write the undulator length as L„; = N„,Ajy, where N„; is the number of 
undulator periods. With the help of Eq. (8) we obtain an expression, valid 
in the far zone: 



U,/2 



exp [zOr] exp 



.codh 

I 

2c 



— sin (k^z') ex + 6 
7 



Here 



(13) 



CO 



[ley 



+ 



2c 



z' - 



Kdx CO 
y k^jC 



cos{k^„z') - 



872 k^c 



sm{2ky„z') . 



(14) 



The choice of the integration limits in Eq. (13) implies that the reference 
system has its origin in the center of the undulator as in Fig. 1. As for Eq. (7) 
and Eq. (8), one recognizes two terms in Eq. (13). The first in sin(^:a,z') can 

be traced back to the first (current) term in Eq. (7), while the second, in 6, 
corresponds to the second (gradient) term in Eq. (7). 

Eq. (13) has been derived under the paraxial approximation, that can always 
be applied for ultrarelativistic systems characterized by a large parameter 
7^ » 1. 
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Fig. 1. Geometry for undulator radiation. 
2.3 Resonant approximation in the far zone 

Usually, it does not make sense to calculate the intensity distribution from 
Eq. (13) alone, without extra-terms (both interfering and not) from the other 
parts of the electron trajectory. This means that one should have complete 
information about the electron trajectory and calculate extra-terms to be 
added to Eq. (13) in order to have the total field from a given setup. Yet, 
we can find particular situations for which the contribution from Eq. (13) is 
dominant with respect to others. In this case Eq. (13), alone, has independent 
physical meaning. 

One of these situations is when the resonance approximation is valid. This 
approximation does not replace the paraxial one, based on » 1, but it is 
used together with it. It takes advantage of another parameter that is usually 
large, i.e. the number of undulator periods Nk, » 1. In this case, the integral 
in dz' in Eq. (13) exhibits simplifications, independently of the frequency of 
interest due to the long integration range with respect to the scale of the 
undulator period. 

In the particular case when the frequency of interest is near the fundamental 
resonance frequency 



cOr = Ikujcyl , (15) 

or other harmonics odd multiples of cOr, extra-simplifications can be ex- 
ploited, allowing one to neglect the gradient term in in Eq. (13), as well as 
the constrained particle motion in the Green's function, i.e. the second term 
in cos(A:j(,z') in the phase Eq. (14), that corresponds to the term in 6 • r'i_{z') 
in Eq. (8). This leads to horizontally polarized radiation and to azimuthal 
symmetry of the field. It should be stressed that odd harmonics constitute a 
particular case. Neglecting the gradient term and the constrained particle's 
motion in the Green's function does not coincide with the application of the 
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resonance approximation, understood as exploitation of the large parame- 
ter Nu, » 1. For example, resonance approximation can be used to study 
even harmonics but, as shown in [23], in that case the gradient term and the 
constrained particle's motion in the Green's function must be retained. 

In this paper we will be interested in frequencies near the first (funda- 
mental) harmonic cOr, so that the above mentioned extra-simplifications can 
be exploited. Let us show how this can be done, and let us discuss how 
radiation characteristics (polarization and symmetry) are related to the pos- 
sibility of neglecting gradient term and constrained particle's motion in the 
Green's function. First, we can specify "how near" co is to co,. by introducing 
a detuning parameter C, defined as 



C 



CO 



lylc 



-k -—k 



(Or 



(16) 



Here co = cOr + Aco. Eq. (14) can thus be written as 



/, ^ co0^\ , Kdr CO ,. CO . ,. ,^ „. 

Or = /c,„ + C + -— z' £ cos(fc,„z') - — — sm(2fc,„z') , (17) 

so that, in all generality, the field in Eq. (13) can be written as 



E_L = exp 
X exp 



.CO eh 



Ic 



Uvl2 

^ ^ ciz'|^[exp(2zl,„z')-l]e, + 0exp(zl,„z')j 



KO, CO 

y h„c 



CO 

cos{k^z') - — sin(2fc^z') 



(18) 



As first proposed in [25|1 one may use the Anger-Jacobi expansion: 

oo 

exp [ia sin{ip)] = ^ Jp{a)exp[ipip] , (19) 



where /p(-) indicates the Bessel function of the first kind of order p, to write 
the integral in Eq. (18) in a different way: 



Ej_ = exp 



.CO eh 



3^ Yj Jm{u)}niv)exp 



inn 
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X ^ dz' exp 

-UI2. 



i\C + 



Ic 



^ [exp {2ik^z') - 1] + exp HKz') ^ 



X exp [z(n + 2m)kwz'\ , 
where 



(20) 



(21) 



Up to now we just re-wrote Eq. (13) in a different way. Eq. (13) and Eq. (20) 
are equivalent. Of course, definition of C in Eq. (16) is suited to investigate 
frequencies around the fundamental harmonic but no approximation is 
taken besides the paraxial approximation. 

Whenever 

C + — « fc^ , (22) 



the first phase term in z' under the integral sign in Eq. (20) is varying slowly 
on the scale of the undulator period Aa,. As a result, simplifications arise 
when Na, » 1, because fast oscillating terms in powers of exp[iky,z'] ef- 
fectively average to zero. When these simplifications are taken, resonance 
approximation is applied, in the sense that one exploits the large parameter 
N-u, » 1. This is possible under condition (22). Note that (22) restricts the 
range of frequencies for positive values of C independently of the obser- 
vation angle 6, but for any value C < (i.e. for wavelengths longer than 
Ar - c/cor) there is always some range of 6 such that Eq. (22) can be ap- 
plied. Altogether, application of the resonance approximation is possible 
for frequencies around cOr and lower than coj.. Once any frequency is fixed, 
(22) poses constraints on the observation region where the resonance ap- 
proximation applies. Similar reasonings can be done for frequencies around 
higher harmonics with a more convenient definition of the detuning param- 
eter C. 

Within the resonance approximation we further select frequencies such that 



■ — - «; 1 , i.e. |C| «: . (23) 

Note that this condition on frequencies automatically selects observation 
angles of interest 6^ «: l/yl- In fact, if one considers observation angles 
outside the range 6^ ^/ylf condition (22) is not fulfilled, and the inte- 
grand in Eq. (20) exhibits fast oscillations on the integration scale L^;. As a 
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result, one obtains zero transverse field, E± = 0, with accuracy 1/Na,. Under 
the constraint imposed by (23), independently of the value of K and for 
observation angles of interest 6^ ly\, we have 

, , K\ey\ w Aa;\ 2^f2K ,^ , ^ ,^ , 

\v\ = —7— = 1 + — , fz Ox < fz OA « 1 . 24 



This means that, independently of K, \v\ <^c 1 and we may expand ]n{'^) in 
Eq. (20) according to /„(u) ^ [2""/r(l + n)] u", !(•) being the Euler gamma 
function 



r(z) = J ^'"^ exp[-f] . (25) 



Similar reasonings can be done for frequencies around higher harmonics 
with a different definition of the detuning parameter C. However, around 
odd harmonics, the before-mentioned expansion, together with the appli- 
cation of the resonance approximation for Ny^ » 1 (fast oscillating terms in 
powers of exp [zl„,z' ] effectively average to zero), yields extra-simplifications. 

Here we are dealing specifically with the first harmonic. Therefore, these 
extra-simplifications apply. They enforce a stronger version of the resonance 
approximation allowing one to neglect both the constrained motion in the 
Green's function and the gradient term in the expression for the field (re- 
spectively, the term in cos(fca?z') in the phase of Eq. (18) and the term in 
in Eq. (18)). First, non-negligible terms in the expansion of ]n{v) are those 
for small values of n, since /„(y) ~ v", with \v\ ^ 1. The value n - gives 
a non-negligible contribution Jo{v) ~ 1. Then, since the integration in dz' is 
performed over a large number of undulator periods Na, » 1, all terms of 
the expansion in Eq. (20) but those for m = -1 and m = average to zero 
due to resonance approximation. Note that surviving contributions are pro- 
portional to K/y, and can be traced back to the current term in only, while 

the gradient term in in Eq. (20) averages to zero for n — 0. Values n — ±1 
already give negligible contributions. In fact, }±i{v) ~ v. Then, the term in 
Cx in Eq. (20) is v times the term with n = and is immediately negligible, 
regardless of the values of m. The gradient term in d would survive aver- 
aging when n = 1, m = -1 and when n = -1, m = 0. However, it scales as 
-» 

dv. Now, using condition (23) we see that, for observation angles of interest 
02 «: l/yl, |0| ~ ( V2 K / V2 + K2 ) y^02 ^ Therefore, the gradient 
term is negligible with respect to the current term in ex for n = 0, that scales 
as K/y. All terms corresponding to larger values of \n\ are negligible. 
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Summing up, all terms of the expansion in Eq. (19) but those f orn = and 
m = -lorm-0 give negligible contribution. After definition of 



Aj 



coK^ 



coK^ 



that can be calculated at co = co^ since |C| «; ku,, we have 



(26) 



^ KcoeAjj 



.0)6^2 



exp 



-1.^/2 



r 



(27) 



yielding the well-known free-space field distribution: 



.codh 



sine 



^(c + 



2c 



(28) 



where sinc(-) = sin(-)/(-)- Therefore, the field is horizontally polarized and 
azimuthal symmetric. It should be clear that result in Eq. (28) is valid for 
arbitrary values of K, that may assume values much larger than unity, as in 
the long-wavelength range of the FIR undulator line at FLASH. 

Condition (22) defines the applicability region of the resonance approxima- 
tion in free-space in terms of the detuning parameter C and the observation 

angle 6. Condition (23) defines the range of the detuning parameter C where 
a stronger version of the resonance approximation applies for odd harmon- 
ics. It imposes a stronger constraint on the detuning parameter C but it 
eliminates constraints on the observation angles because, as said before, 
observation angles outside the region of interest 0^ «; l/yl simply return 
zero field. Thus, Eq. (28) can be used without restrictions on 0. Sometimes 
we will say that (23) enforces the "strong" resonance approximation, when 
the extra-simplifications we just discussed apply. 

The constrained motion term in the Green's function and the gradient term 
in the equation for the field are strictly related to polarization and azimuthal 
symmetry of the field. When the constrained motion in the Green's function 
can be neglected, also the gradient term can be neglected and the field 
is horizontally polarized and azimuthal symmetric. There are situations 
when this cannot be done. For example, when studying the long-wavelength 
region when (23) is not valid, the field looses these properties. However, if 
(22) holds, resonance approximation can still be applied, in the sense that the 
large parameter » 1 can be exploited to simplify the integration in dz' . 
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In the case of even harmonics, that has been studied in reference |23l] , one 
can use restrictions (23) and expand the Bessel function similarly to what 
has been done here. However, in this case, both constrained motion term in 
the Green's function and gradient term in the equation for the field are not 
negligible. Therefore, the field is not horizontally polarized, nor azimuthal 
symmetric. 



2.4 Undulator radiation in the near zone 



As has been described in [16], SR fields can be treated in terms of Fourier 
optics. In particular, radiation from an ultra-relativistic particle in an un- 
dulator can be interpreted as radiation from a virtual source placed in the 
undulator center, i.e. at z = 0, which produces a laser-like beam. The con- 
cept of virtual source has already been proposed (see | |27tl , and (29']). 
In fld\ we developed this concept up to an analytical description of the 
virtual source that is similar, in many aspects, to the waist of a laser beam. 
In the case of undulator radiation in free-space it exhibits a plane wave- 
front. It is completely specified, for any given polarization component, by 
a real-valued field amplitude distribution and can be determined from the 
knowledge of the far-zone distribution in Eq. (28) through the relation (see 
(kl): 



£x(0,r-l) = — j ddexp 



ia)\d\^ 
2c 



Ej_{d) exp 



—r^ ■ e 
c 



(29) 



where the integration in dO is understood to be performed over the entire 
plane spanned by 6. The integral in Eq. (29) can be calculated analytically 
at perfect resonance (C = 0) yielding the field distribution at the source: 



£x(o,fj = - 



iAjjedgCO 



2c2 



7i-2Si 



'rjcv^ 



(30) 



where Si(-) indicates the sin integral function 



Si(z) 



z 

-I 



dt 



sm{t) 



(31) 



The quantity in Eq. (30) is a scalar because, as we have seen, the field is 
horizontally polarized. Numerical investigations are possible off-resonance. 
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for C 0. Once the field at the virtual source is known, the field at other 
longitudinal positions, both in the far and in the near zone, up to distances 
to the sources comparable with the radiation wavelength, can be obtained 
with the help of the Fresnel propagation formula: 



£±(z,rL) = 



l(x) 

Inez 



J dr\ £x(0,r\ 



)exp 



10) r I 



2cz 



(32) 



yielding, at any position z 

KcoeAjj 



E± (z, rx) = 



2c^y 



Ei 



, 2zc LmC I 



Ei 



> 2zc + LjnC J 



(33) 



Here Ei(-) indicates the exponential integral function 



Ei(z) Jdt f^^il , (34) 

-z 

where the principal value of the integral is taken. The far-zone limit of this 
expression is reached for z » L„; and reads 



~ K(veL-a,Ajj 
£±(z, 0) = exp 



lyc'^z 

that coincides with Eq. (28) at C = 0. 











' Ic 


sine 


4c 





(35) 



3 Paraxial Green's function with boundary conditions 



In the previous Section 2 we considered the problem of characterizing radi- 
ation from an ultrarelativistic electron moving on an arbitrary trajectory in 
free-space using Paraxial Maxwell's equations. 

We now wish to generalize that treatment, to describe radiation from an 

ultrarelativistic electron moving on a given trajectory inside a metallic pipe 
of arbitrary cross-section. An example is depicted in Fig. 2, where n is a 
vector field defined on the boundary surface S, such that \n\ = 1. At any 
point, ft is orthogonal to S and points outwards. 
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metallic boundary surface 

s 



Waveguide 
cross-section 




Electron moving 
inside the pipe 



Fig. 2. Geometry of the problem. The electron is moving inside the pipe. 

As we have discussed at the end of Section 2.1 by expanding the far-zone 
field distribution in free-space in terms of plane waves, the paraxial approx- 
imation applies when the transverse wave number k_i_ obeys k]_c^/co^ ^ 1. 
This is no restriction for ultra-relativistic systems with 1 /y ^ «; 1 . In fact, non- 
negligible components of the angular spectrum are those for k]^c^/co^ «: 1. 
We can take advantage of this relation also in case a waveguide of typical 
geometrical dimension R is present (for a circular waveguide R indicates 
the radius). The difference with respect to the free-space case is in a discrete 
number of modes, corresponding to a discrete range of values that k± may 
assume. The concept of radiation formation length must be ascribed to each 
propagating mode, but for ultra-relativistic electrons Lf ^ A holds for each 
mode. One concludes that paraxial approximation holds when k^c^/co^ <^ 1 
for each mode. The same conclusion can be drawn by remembering that the 
paraxial approximation is based on the fact that the electron velocity must 
be close to the phase velocity of light. In the waveguide, the phase velocity 
of light is larger than c, i.e. t^ph = c(l -c^k]^ /cu^)"^^^. As a result, in order for the 
electron velocity to be near to Vph one must require that the electron velocity 
is close to c and additionally that k^c^lap- <^ 1. This means that the paraxial 
approximation is applicable for frequencies well above the cutoff frequency 
of a given mode, co » cucutoff - ck^. When co » tUoitoff the waveguide is said 
to be overmoded. Thus, the paraxial approximation is applicable when the 
waveguide is overmoded. Since the smallest cutoff frequency (that for the 
fundamental mode) is about cucutoff ~ c/R, it also follows that the waveg- 
uide is overmoded for A «; R. We assume that the last condition is verified 
through this paper. 

Since A <^R the paraxial approximation can be applied and is convenient to 
introduce, as for the free-space case in Section 2, the slowly varying envelope 

of the Fourier transform of the transverse electric field, "the field" Ex(z, r^_). 
Here (z, r^_) indicates a point inside the waveguide, as shown in Fig. 2. The 
paraxial equation for the field is identical to that for the free-space problem. 
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Eq. (3), where the differential operator T> is defined in Eq. (4), and /(z, r^) in 
Eq. (5). 

We only need to account for the presence of the vacuum pipe by specifying 
proper boundary conditions. In the case of perfect metal (we will treat wall 
resistance effects in Section 6), the electric field must be orthogonal to the 
boundary surface S, i.e. 

(n^Ey-HyEX =[nxE_L)| =0 (36) 



and 



(37) 



Maxwell's equation for the curl of the magnetic field implies that Eq. (37) is 
equivalent to 



VxH 



ZJ 



1=0. (38) 

Is 



Here H, analogously to E, is the slowly varying envelope of the magnetic 
field in the space-frequency domain, and will be referred to as the "magnetic 
field" . Since paraxial approximation applies, the transverse components of 
the magnetic field and Hy are related to the electric field by E^ - Hy and 
Ey - -Hx so that Eq. (38) finally yields 



(dE, dEy^ 
dx dy 



Vx-E^i =0 (39) 



Eq. (36) and Eq. (39) are boundary conditions in the paraxial approximation. 
To the best of our knowledge, they are only discussed in |j2l|]. At first glance 
the second condition, Eq. (39), is not obvious. In paraxial approximation Eq. 
(37) looks automatically satisfied, so that one might be led to use Eq. (36) 
only. Actually, Eq. (36) poses conditions on the field amplitude while Eq. 
(39) poses conditions on the field derivatives, so they are both necessary. 
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Summing up, we should jfind solution to the problem 



£±(z,rx) 



(nx£x)| = 
(Vx-£x)| =0, 







(40) 



where the operator D and the vector function / are defined in Eq. (4) and 
Eq. (5) respectively. We proceed by using a Laplace transform technique, 
that allows one to dispose of the partial derivative with respect to z in 1) in 
the first equation in (40). 

First we define Laplace X and inverse Laplace transform of a function 
g(z,rx)as 



00 

g(p, r^) = £ [g{z, rl)] ip) = J g{z, rl) exp [-pz] dz 



(41) 



with Re[p] > 0, and 

g{z, fj_) = [gip, fx)] (z) = J g(p, r±) exp [pz] dp , 



6+100 



(42) 



6-ioo 



where 6 is a real number larger than all the real parts of the singularities of 

giv)- 

Consistently applying a Laplace transformation to the equation set (40) we 
obtain a redefinition of the problem in terms of the Laplace transform of the 



field E_L(p,r±) = £ 



E±{z,r_L) 



ip): 







(Vx-£x)| =0, 



= {2ico/c)EA0,n)+f{p,ri.) 



(43) 



where / is the Laplace transform of / and 
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(44) 



Note that the presence of the initial condition Ej^(0, rj_) in Eq. (43) refers 

to the possibility of introducing an external field into the system. In what 
-» 

follows we set Ej_{0, r±) = 0. 

Now, suppose that we find a tensor of components G" such that 



£« = j Q;{A,r\,p)f{r\,p) dr' , 



(45) 



where we indicated tensor and vector components with Greek indexefl In 
this case, the inverse Laplace transform of G^, that will be written as G"^, 
is the tensor Green's function for the problem (40), inclusive of the proper 
boundary conditions. As a resu 10 



z 

'(rl, z) = J dz' J dr\ Gl (rl, r\, z - z') ff" {r\, z') , 



(46) 



where, as in Section 2.1, we integrate up to z' = z because the radiation 
formation length for z - z' < is very short with respect to the case z - z' > 

0. Summing up, we first have to find G, then to apply a Laplace inverse 
transform in order to get G and, finally, to solve for E". 

We start by specifying the eigenvalue problem associated with the problem 
set (43), that is 

'5[f,(r-l)] = A/,(rl) 

< X ^/)|^ = (47) 
(V,.f,)j =0. 

V Is 



Here the index / is understood to identify any possible vector solving the 
equations set in (47). Posing Ay = licop/c-Aj yields an equivalent eigenvalue 

^ Note that high and low indexes are interchangeable, as the metric tensor is 
represented by the identity matrix. 

-* 

Formally we are translating the initial point z = 0, where E_l(0, r^) = 0, at z — > 
— oo. 
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problem involving Helmholtz's equation: 

'vif/rl) + A/y(r-l) = 

.(«x^i)|^=0 (48) 
(Vx-F,)| =0. 

Is 



It can be verified that boundary conditions are homogeneous, so that the 
domain of the Laplacian operator is, in our case, the vector space of twice 
differentiable (square integrable) functions obeying boundary conditions 
in (48). The Laplacian operator defined in this way can be shown to be 
self -adjoint with respect to the inner product defined as 

{f'i)^ j r-gdr\- (49) 

s 



This means that for any chosen / and g in the domain of the Laplacian 
operator we have 

(/.Vi#y = (^-,Vi/), (50) 

where asterisks indicate complex conjugation. As a result, eigenvalues Aj are 

real. Eigenvectors Fy referring to different eigenvalues are orthogonal and 
span over the entire domain. The spectrum is discrete, and eigenvectors can 
be chosen real. We normalize the eigenvectors Fj so that 



(F;vFy} = l. (51) 
As a result we can write the decomposition 

^ L 2icop/c-Aj ' ^^^^ 



that is analogous to the approach for solution of a Sturm-Liouville problem. 
Since eigenvectors can be chosen real we will do so, and from now on we 
will avoid complex conjugation everywhere. 

Helmholtz's theorem (also named the fundamental theorem of vector calcu- 
lus) states that any vector field (in our case, the two-dimensional vector field 



24 



Fj) that is twice continuously differentiable and vanishing rapidly enough 

at infinity, can be split in the sum of two vector fields; the first, Vj, irrota- 

-* -* -> -> 

tional, i.e. V x V, = and the second, V^, divergenceless, i.e. V • = 0. 

Considering two independent scalar (and real) functions i/'J^ and we 
take Vi = d^ip™e, + dyip™ey and = dyip'^^e, - We thus obtain the 

following representation of Fj in terms of scalar functions ip"^^ and ip™: 



Fj = [el. + ^.^™) + 4 (-d.^pf + <9,^™)] . (53) 

In the following it will be clear that \p^^ and \p™ are physically related, 
respectively, to transverse electric and transverse magnetic field modes. At 
the present time, however, these functions are only introduced for mathe- 
matical convenience. In fact, by substitution in (48), the vectorial eigenvalue 
-> 

problem for Fj can be reformulated in terms of the following scalar problem 
\lxP™{r^) + Aj''™ »/'™(rx) = 

< "^•(V.^f)|^=0 (54) 

(0| =0- 

V Is 



Boundary conditions in Eq. (54) do not allow for hybrid TE-TM modes, as 
in the case of walls with finite resistance (see Section 7). In other words, 
solutions ip^^ and ip™ are not degenerate. This can be seen by inspection of 
Eq.(540 ' 

It follows that the particular choice of boundary conditions allows a choice 
of eigenvectors Fj different from that in Eq. (53), that gives the same results 
but allows separation of TE and TM modes: 



. / ^ " . (55) 

^ In fact, suppose there exist two non-zero solutions if^^ and ip™ sharing a com- 
mon eigenvalue A. One would have V^(i/'J^ - ip™) - 0. However, on S, ip™ — 
and thus, on S, V^i/'J^ = 0. For this to be in agreement with Helmholtz's equation 
in (54), one must have A = 0, hence V^i/^J^ = and V'^ip™ - everywhere inside 
the waveguide. Then, since ip™ = on S it also must be zero everywhere, and 
degeneracy is impossible. 
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We formally adopt this new eigenvector set, because this is compatible with 

rill" T^l jf" 

the following normalization condition for ifjj ' : 



- TEJM 



' = 1 . (56) 



The reason why the eigenvector set (55) and the normalization condition Eq. 
(56) are desirable is only formal: this choice yields traditional results known 
in literature for the eigenvalue problem (54) when a particular geometry of 
the waveguide is fixed, as in the next Section 4. 

Solution of the problem (54) depends on the setup geometry and cannot be 
further specified unless the surface S is given. Later on we will consider the 
case of a circular waveguide. At this time we assume that the problem is 
solved in some particular geometry, and that if^.^'™ are fixed. As a result Fj 

are known through Eq. (53) and the tensor G"f^ is completely specified. In 
order to find the tensor Green's function, we should find the inverse Laplace 
transform of G"^, that is 

6+ioo 

G«^^(z,fl)^X-i[G«^(p,rl)](z)= j G"^p,i'^)exp[pz]dp. (57) 

6-ioo 



Here 6 is a real number positioned, on the complex plane, on the right-hand 
side of all poles pj . These are located at 

Pi = -i^, (58) 



as it can be seen by inspecting Eq. (52). 

From Eq. (52), and from the fact that Fj - Fj{r±), we see that the dependence 

of G on p is only in the denominator of the right hand side of Eq. (52). As a 

result, G — > uniformly as p — > oo. It follows that Jordan's lemma applies 
and we obtain 



G«^ = 2],{Res[Gexp(pyz)]f^ (59) 
where 
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(Res [g exp[pyz]]) 



lim 

-icAjl{la)) 



icA 



|p + _i|G«^(p)exp[p,z]j 



2ico 



Ico , 

F" (rl)F^(r\)exp 



2a; 



(60) 



Substituting Eq. (53) in Eq. (60), and Eq. (60) in Eq. (59) we obtain the 
following explicit expression for the tensor Green's function 



icAf 



(rJ-4^.V^f (rl)] 



+ 



—T 



exp 



icA™ 



[.-;^,^™(r-l) + 4<9,V^™(rl)] 



(61) 



where notation indicates the tensor product. It should be noted that the 
result in Eq. (61) coincides with that obtained in [210. Once the tensor 
Green's function G is known, the slowly varying envelope of the electric 
field is given by 



z 

j dz' {^t/Az')G; {A,r\{z'),z- z') + d'fil (r-l,r\(z'),z - z')j 

2' 

ico r 



dz 



where we made explicit use of Eq. (2), and derivatives d'^ are taken with 
-* 

respect to r'_L . 



^ Numerical factor and sign differences are due to different definitions for the 
Green's function used in 12111 compared to the present paper. In that reference a 
Green's function is first calculated starting with the rigorous solution for the eigen- 
functions of the waveguide, without applying paraxial approximation. Only in 
a second step such approximation is applied. In this paper instead, we directly 
derive the paraxial Green's function for an overmoded waveguide using approx- 
imate boundary conditions on the waveguide walls (obtained by application of 
the paraxial approximation) and a Laplace-transform technique. We thus directly 
solve the initial value problem in paraxial approximation. The present derivation 
is shorter and more straightforward than that in |,21. ]. 
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4 Green's function for a circular waveguide 



We now restrict our analysis to a particular geometry of the waveguide. 
Namely, we select a circular cross-section of radius R. In order to explicitly 
compute the tensor Green's function in Eq. (61) we need to solve the problem 
in (54) with this particular choice of geometry. 

For both TE and TM modes, we need to solve Hehnholtz's equation. Be- 
cause of the particular symmetry in the case of a circular waveguide, it is 
convenient to use transverse polar coordinates so that rx is identified by 
radial distance r - ^x^ + and angle (p - arctan(y/x). The equation set in 
(54) yields 



dr^ r dr d(\P- 



,TE,TM 



TE,TM ,/,TE,TM _ 



= 0, 



(63) 



with boundary conditions 



TE\ 



dr 



= and 



lr=R 



= 0. 



(64) 



Solutions of Eq. (63) are 



/ TE,TM ^ 
'mkl 
, TE,TM 
V 'mlO. ) 



sin(m(^) 

^ cos{m(p) ^ 



(65) 



where the eigenvalues A^f '™ and the normalization coefficients A^™ are 

o mk mk 

specified by the TE and TM boundary conditions in (54). 



Switching to notation y.^]^ = -^A^ R and v^k = -^A™ R one defines j^mk and 
Vmk respectively as the roots of equations 



J'm (^mfc) = 

and equations 



(66) 



Jm {Vmk) = . (67) 

The correspondent normalization coefficients for the TE and for the TM 
modes are given by 
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and 



'"^ y n V„k Jm-l{Vmk) 

with ao = 1 and = 2 for m > 1. 

This means that jfixing an eigenfunction Fy corresponds to the specijfication 
of two positive integers m > and k > 0, and either a sin(-) or a cos(-) 
solution. 

We should now substitute Eq. (65) in the expression for the Green's function, 
Eq. (61). Before doing so, however, we need to express partial derivatives of 
functions i/^^. ' with respect to x and y in terms of partial derivatives with 

respect to r = ^x^ + and (p = arctan(y/x). This can be done with the help 
of relations 

d d \ d 

- = cos{cp)--sMcP)-^ (70) 



and 



5^=smW5^+cosW-5;^. (71) 



Using Eq. (70) and Eq. (71) and taking advantage of: 

(72) 

and 



= \ [Jm-l (a - Jm.l m (73) 

we obtain, from Eq. (61), the following expression for the tensor Green's 
function in the case of a circular waveguide: 



29 



oo oo 



ic{z - z') 2 



X < 



+ 



\2R/ 

m=0 fc=l 

]rn-i{^mkrlR)cos{{m - + /„,+i(fi„^r/R)cos[(m + 
-}m-iil^mkr/K) sin[(w - 1)(/)] + Jm+iil^mkr/R) sin[(OT + 
Jm-iifimkr'/R) cos[(m - 1)(^'] + ]m+x{[imkr' IR) cos[(m + 1)(^'] 
-/m-i(iUm)t?'7-R)sin[(m - + /^+i(|U^fcr7R)sin[(m + 

-Jm-ii^hnkr/R) sin[(m - l)cj)] - }rn+i{^mkr/R) sin[(m + l)(p] 

-]m-MmkrlR) COS[(m - 1)(/)] + Jm+lil^mkr/R) COs[(ot + 1)(/)] 

-Jm-ii^mkr'/R) sin[(m - - Jm+i{^mkr' /R) sin[(m + l)^)'] 
-Jm-iifimkr'/R) cos[(m - 1)(^'] + J^+iijimkr'/R) cos[(m + 



oo c« 



ZC(2 - Z') 2 



X < 



+ 



m=0 
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Jm-i{Vmkr' /R) sin[(m - 1)(^'] - }m+i{Vmkr' /R) sm[{m + l)(p'] 

Jm-i{vmkr'/R)cos[{m - l)(p'] + }„,+i{vrnkr'/R)cos[{m + l)(p'] 

}m-i{Vmkr/R) cos[(m - l)(p] - Jm+i{Vmkr/R) cos[(m + l)(p] 
-Jm-i{Vmkr/R)sm[im - l)(p] - }m+i{v,nkr/R)sm[{m + l)(p] 

Jm-i{Vmkr' /R) cos[(m - l)(f)'] - }m+i{Vmkr' /R) cos[(m + 
-Jm-i{vmkr' /R) sm[{m - l)(p'] - frr,+i{vmkr' /R)sm[{m + l)(p'] 



.(74) 



In order to verify the correctness of Eq. (74) we study the free-space limit, that 
corresponds to the limit for large values of the waveguide radius R — > oo. 

Since we are interested in characterizing the fields over a finite transverse 
direction, and since sources have a finite transverse size, the limit R — > 
oo allows one to substitute Bessel functions in Eq. (68) and Eq. (69) with 
asymptotic expressions for » 1. First remember that 



w^)«y^cos(e - — --], 

One sees that Eq. (66) and Eq. (67) are satisfied by 



(75) 
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(tn 1 \ (tn 3 \ 

y + + - I , Vmfc = 71 ( — + + - j , (76) 

yielding 



/mW--^, /Ll(^-fc)--^- (77) 

As a result we obtain 



(78) 



Substituting Eq. (78) in Eq. (74), posing E, = nk/R and replacing the sum 
over k in Eq. (74) with an integral over d£, we obtain, after cumbersome but 
straightforward calculations: 
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sin[m((/) - (p')] cos[m{(f) - (p')] ^ 
Eq. (79) can also be written as 
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As before, z > z'. The integrals in d£, can be performed yielding 
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xexp 



2c(z - z') 
Using the Anger-Jacobi expansion 



cos[m{(p - (p')] 

cos[m{(p - (p')] ^ 



(81) 



/o(C) + 2j^ i-^^UQ cos{m{cp - (/)')] = exp[-zC cos((j^ - (/)')] 



(82) 



m=l 



we find that in free-space the tensor Green's function Eq. (74) reduces, for 
both vertical and horizontal polarization components, to the scalar Green's 
function in Eq. (6): 



^ 47i(z - z') 
as it must be. 



ICO- 



r_L - r'_L Y 
2c(z - z') 



(83) 



5 Radiation of wiggled electron in circular waveguide 



We now want to apply Eq. (74) and Eq. (62) to study the case of an undulator 
in presence of a circular waveguide in the case of a planar undulator. As 
we have seen before, the transverse electron velocity is taken along the hor- 
izontal axis and is specified in Eq. (9). The longitudinal Lorentz factor yz{z) 
is not constant, due to oscillatory motion of electrons along the longitudinal 
axis. As a result, Eq. (11) holds. One must substitute Eq. (9) and Eq. (11) in 
Eq. (62) to calculate the field components with the help of Eq. (74). 

Eq. (62) has been derived under the paraxial approximation. In Section 2 
we discussed the region of applicability of the paraxial approximation in 
free-space, concluding that it can always be applied for ultra-relativistic 
systems when l/y^ «: 1. In Section 3 we extended our consideration to 
the case when a waveguide is present. We concluded that, for an ultra- 
relativistic system, the paraxial approximation holds when the waveguide 
is overmoded, i.e. R » A. In the case of radiation of wiggled electrons this 
condition is necessary but not sufficient. This can be seen straightforwardly 
for our practical case of interest K > 1, with the help of a geometrical 
argument. From Eq. (9) it may be seen that electrons wiggle with a wiggling 
amplitude 
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2KAy, 



(84) 
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We require R> > f^A » A. In fact, when R » A we can still discuss in 
general about an overmoded waveguide and, thus, about paraxial approx- 
imation, but in order to discuss about radiation from a wiggled electron we 
must impose that the waveguide radius is larger than the amplitude of the 
wiggling motion. It is interesting to note that a different reasoning holds for 
K <^ 1, leading to a similar conclusion. Namely, for K «: 1 we can define 
an inertial frame moving with the same average longitudinal velocity of the 
electron, characterized by the average longitudinal Lorentz factor fz. In this 
frame, the wiggling motion is non-relativistical and the undulator field is 
similar to that of an incoming plane wave with wavelength Ayjiy. We are 
actually dealing with a Thompson-scattering problem, where the electron 
radiates as a dipole at frequency a* = cy/Aa,. Such frequency must be higher 
than the cutoff frequency cucutoff ~ c/RsothatR > A^^/y. For frequencies near 
resonance this requirement translates to R > yA. As a result we may state, 
independently of K, that it makes sense to deal with the case R > yzA » A. 

As has been seen in Section 2 in the case of undulator radiation in free-space, 
the resonance approximation constitutes a further approximation that can 
be applied together with the paraxial one. We concluded that resonance 
approximation allows simplifications by exploitation of the large parameter 
Na, » 1. We also discussed how, for odd harmonics and when condition (23) 
holds, further simplifications take place, allowing one to neglect gradient 
term in the field equation and constrained motion in the Green's function, 
giving a stronger type of resonance approximation. This strong resonance 
approximation can be applied not only in free-space but also in the presence 
of a waveguide, exactly as the paraxial approximation can. In the case of a 
circular waveguide with radius R Eq. (74) must be taken as Green's function 
instead of the free-space Green's function, Eq. (83). In the limit for » 1 
and when certain constraints on the parameter space are verified, similarly 
as for the free-space case, one may neglect the gradient term in Eq. (62), as 
well as the constrained particle motion in the Green's function Eq. (74). 

Let us discuss these simplifications in detail. We may use Anger-Jacobi 
expansion to obtain from Eq. (62) a representation of the field that is the 
analogous of Eq. (20) in free-space. Namely, for any value k, the current 
terms in the integrand in dz' in Eq. (62) are proportional to 
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(85) 



while the gradient terms are proportional to 



f g(z') = ^ exp[i{2p + l)ku,z'] exp 




(86) 



where C = f^mit or C = i^mic/ <5/ can be either cos[l(/)'(z')] or sin[l(/)'(z')] arid 
parameter w has been defined in Eq. (21). As one may see by inspection, the 
role of the index p in Eq. (85) and Eq. (86) is the same as the role of m in Eq. 
(20), while the role of the index m in Eq. (85) and Eq. (86) is similar to the 
role of n in Eq. (20). 

Comparison between Eq. (85) or Eq. (86) and the integrand in Eq. (20) shows 
a striking similarity between the role of the observation angle 6 in free-space 
and, for each mode, of the quantity cC/icoR). This relation should not be 
surprising. As we have seen at the beginning of Section 2.1, after expanding 
the far-zone field distribution in free-space in terms of plane waves one 
sees that the transverse wave numbers of interest k± obey k]_c^/oj^ <c 1, i.e. 
the paraxial approximation can always be applied. In fact, the plane wave 
components of the angular spectrum for which paraxial approximation is 
applicable, are also the only relevant components of the spectrum, and are 
those that can be seen at observation angles 6^ <«c 1. The quantity cC,/{a)R) 
is the analogous of the free-space propagation angle of a spacial Fourier 

component, kj_c/co, that is allowed to vary continuously across the reciprocal 
space. The difference with respect to the free-space case is in the presence 
of a discrete number of modes, corresponding to a discrete range of values 
that cC,/{coR) may assume. All we have to do to show the applicability of the 
resonance approximation is to rely on the analogy between 6 and cC,/{a)R), 
and to follow the same derivation for the free-space case. 

First, as already said, Eq. (85) and Eq. (86) are analogous to the integrand of 
Eq. (20). 

Second, condition (22) can be stated, for each mode, as 

C + ^«/c., (87) 



that allows to exploit the large parameter » 1 by averaging powers of 
exp[ffca,z'] over a large number of undulator periods. 
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Third, condition (23) allowed to neglect constrained motion in the Green's 
function and gradient term in the field equation in free-space. This enforced 
the strong resonance approximation around the fundamental (and other 
odd) harmonic. In the presence of a waveguide condition (23) remains un- 
varied: 



« 1 , i.e. |C| « . (88) 



COr 



Note that this condition on frequencies automatically selects propagation 
angles of interest c^C^ / {2coR^) <^l/fl.In fact, one can consider propagation 
angles outside the range c(} I {IcoB?) «: 1 ly\, but in this case Eq. (85) and Eq. 

(86) exhibit fast oscillations on the integration scale L-^ because condition 

(87) is no more fulfilled. As a result, with accuracy l/Nu,, one obtains zero 

transverse field, Ej^ - 0. These remarks are analogous to the free-space case, 
discussed in Section 2.3. Also similarly to the free-space case, the argument 
in the Bessel functions /m±i(-) and /„(•) in Eq. (85) and Eq. (86) is small for 
propagation angles of interest c<^l{lcoF?-) «; l/y: 



,-2 



Cr'(z') Cr. _ 2V2K cC 



This means that, independently of K, we may expand the Bessel functions 
/m±i(0 and ]m{-) in Eq. (85) and (86), exactly as we did in Section 2.3 with 
Eq. (20). As for the free-space case, this expansion can be performed also 
around even harmonics. However, for the first (and odd) harmonic extra- 
simplifications hold. Using the same kind of reasonings as in Section 2.3 
and making explicit use of condition (88) we see that, after integration in 
dz' along a large number of undulator periods N.^, » 1, the only surviving 
contributions are of current type for m = 1 and either /? = -1 or /? = 0. It 
reads: 



F,(z') = ^/p(M)exp[z(c-K^C^)z' 



(90) 



When condition (88) is valid, i.e. |C| ^ k^,, we can neglect, for the first 
harmonic, gradient term and constrained motion in the Green's function, 
exactly as in the free-space case. Note that Eq. (90) is valid also for propa- 
gation angles outside the region of interest cC^/{2coR^) «c l/fl- In this case 
in fact, once the integration in dz' is performed, one obtains a negligible 
contribution to the field, with accuracy 1/Na,. 

Let us consider the propagation direction of interest cC,^/{2ijoR^) <c l/fl- 
It can be seen that, in order to obtain non-zero contributions, we may 
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take C ~ llZl Thus, our condition on the propagation directions of inter- 
est becomes » ifz^)^- This condition should be stated together with 
the previously found constraint > (fz^)^ ^ meaning that parax- 
ial approximation can be applied and that we can discuss about wig- 
gled electron motion. Note that we always deal with the small parame- 
ter (Ayz)^/(ALa,) ~ XjN-u, ^ 1, where KL-^a is the squared radiation diffrac- 
tion size. It follows that we can discuss about both strong {R} < ALf„) 
and weak (R^ » AL^i,) waveguide influence within the required condition 
R^ » (Afz)^ » A^. It is also interesting to remind that in cases of practical 
interest (this is the case for the undulator beamline at FLASH) we have 
K > 1, so that r^u - IKy^. 

Summing up, Eq. (62) gives 



£« = ^ J rf2'exp[^Cz']G«|,„(^,^^„, (91) 

where Ajj is defined in Eq. (26), and G is given in Eq. (74). We obtain 



E,,{r,(p,z) = 



icoedoA 



k=l 



/o + Ji (fti)c^jcos(2(/)) 

Jo (^it^) - /2 (vi)t^)cos(2(/)) 



-I- 



(92) 



and 



Eyir,(p,z) = -!^^ifiIsin(2(/))£ /^{^'(z)/^ (^^i,^) - :?i;;(z)/2 (vi,^)}(93) 

k=l 

where :?lj'(z) and :?l^(z) are given by 

^ From the phase factor in Eq. (62) we see that there is a maximal integration length 
of interest, characteristic of any system, after which the integrand of Eq. (62) exhibits 
oscillatory behavior. This length is the formation length for the system in free space, 
Ljr. Once L J is known, from the phase factor in (z-z') in Eq. (74) follows that there is 
an upper limit to the values of C, giving non-negligible contributions (with a certain 
fixed accuracy) to the field. In fact, above a certain value of C, the exponential 
functions in (z - z') in Eq. (74) show strongly oscillating behavior in the integration 
variable z'. This upper value is estimated as the ratio between the pipe radius R 
and the radiation diffraction size yfXTf, that is R/ ^JALf. We conclude that values 
of interest for C, are up to order R/ yJALf. This also mean that, in calculations, the 
maximal number /Cmax of modes of interest should be taken as as k^ax ^ R/ ^f^Lf■ 
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exp[-iC^z] 



L„,sinc 



7 V-k 



(94) 



and 



qexp[-fqz] 



,2 n 



Lu,smc 



wo 
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(95) 



where 



2a;R2 



^1^ 

2a;R2 



(96) 



Note that each mode propagates according to an exponential law in z, as 
specified in Eq. (94) and Eq. (95). The field can be propagated forwards 
or backwards, changing the sign of z. The field at position z - can be 
interpreted as radiation from a virtual source placed in the undulator center 
and producing a laser-like beam, similarly as in the free-space case. 

As before, it is instructive to study the free-space limit. Proceeding as in the 
previous Section 4, i.e. using asymptotic expressions for J^i^ik) and J^ivik), 
and substituting the summation over k with an integral in d£„ with = nk/R, 
we obtain the following free-space limit for the horizontal field component: 



Ex{zj,C) = — I <Ssmc 



Jo{r^) exp 



i^^cz 
'~2^ 



(97) 



The integral in Eq. (97) can be solved analytically at perfect resonance (C = 0) 
and for z = 0: 



£.(0,r,0) = 



iAjjedgCo 
2c2 



7z-2Si 



(98) 



Position z = is in the center of the undulator. The meaning of Eq. (98) is, 
in fact, that of the field distribution of a virtual source located in the middle 
of the setup discussed in Eq. (30) of Section 2.4. Similar reasoning applied 
to the vertical component gives back zero, that is in agreement with the fact 
that the free-space field from a planar undulator (with vertical magnetic 
field) is horizontally polarized. 

In free-space, under the applicability of the strong resonance approxima- 
tion (condition (23)), one obtains a horizontally polarized field that is also 
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azimuthal symmetric. This is not the case when one is working under the 
applicability of analogous approximation in the presence of a waveguide 
(condition (88)). The presence of a waveguide alters boundary conditions of 
the problem, that are now formulated mixing horizontal and vertical polar- 
ization components (see Eq. (36)), thus destroying horizontal polarization 
and azimuthal symmetry. 

To conclude this Section we mention that in Appendix A we present a 
derivation of the field in the case of a helical undulator. Comparison with 
the planar undulator case, that we just treated, yields immediately 




where E* are the fields from electrons rotating in opposite directions in a 
helical undulator, as defined in Appendix A. Eq. (99) means that the field 
from a planar undulator can be seen as a superposition of fields from two 
helical undulators where electrons follow clockwise and counterclockwise 
trajectories. 



6 Analysis of results (Exemplifications) 



We now analyze our main results, Eq. (92) and Eq. (93). To this purpose, it 
is convenient to introduce the following normalized units: 



I -1—^—\e 

^ \ Ajjcoeds) 
C = CL^ = 2nN^— 

Z 

^ = r 

r 

r = 



n = — 

q = qi^ = ^ • (100) 

We may write Eq. (92) and Eq. (93) in normalized units as 
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(101) 



and 
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(102) 



J?lj"(z) and ^Az) can similarly be expressed as 



mz) = 



sine 



(103) 



and 



^(z) = 



Clexpi-iOz] 



-sine 



Ifc-'O 
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(104) 



Note that the sinc(-) functions in the expressions for is a direct conse- 
quence of our model of the undulator, where the magnetic field is instanta- 
neously switched on and off, at positions z = -1/2 and z = 1/2 respectively. 
In fact, the sinc(-) functions are Fourier transforms of a rectangular function 
with respect to + C or + C, modelling the instantaneous switch on and 
switch off of the undulator field. The presence of high frequency compo- 
nents in the rectangular function implies the presence of contributions with 
high values of k. In its turn, the sum over k in Eq. (74) can be interpreted as 
a superposition of plane waves propagating at angles /ii^/ Vo (or v^t/ VO) 
in units of the diffraction angle ^JXJh^. Thus, higher values of k correspond 
to the introduction of high spatial frequency components in the expressions 
for the field, Eq. (101) and Eq. (102). 

However, we should account for the fact that our theory applies with a finite 
accuracy related to the use of the resonance approximation. When using 
this approximation we neglect contributions to the field with an accuracy 
of order 1/N„;. The instantaneous switch on and switch off is also valid with 
this accuracy This means that, in the analysis of our results, it does not 
make sense to consider high spatial frequency contributions due to abrupt 
switching of the undulator fields on a scale shorter than the undulator period 
Aa,, because these are outside of the accuracy of the paraxial approximation. 
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We may introduce a spatial frequency jfilter in our expression for the field 
by replacing Eq. (91) with 



oo 

~ 2neco9sAjj C , 



(105) 



where the function S(z') introduces some smoothing of the rectangular un- 
dulator profile on a scale of Ajy. In normalized units Eq. (105) reads 



oo 

£« = -27i J dz'exp[iCz']S{z')G%,^.,^^^. 



(106) 



We model S(z') as a constant function along the undulator length with 
exponentially decaying edges on a typical distance A: 



exp [-(z' + 1/2)2/(2A2)] for z' < -1/2 
S(z') = ^ 1 for - 1/2 < z' < 1/2 

exp [-(z' - 1/2)V(2A2)] for z' > 1/2 

It follows that Eq. (103) and Eq. (104) should be replaced by 



(107) 



(108) 



and 



C^expf-zC^zl 

where y{S(z'), (C^'^ + C)} is the Fourier transform of the function S with 
respect to (C[''^ + C). The introduction of the function S(z') introduces a sup- 
pression of higher frequency components of fiS}, and corresponds to a 
suppression of higher spacial frequencies in the field distribution, i.e. to a 
smoother field distribution. For a fixed value of Q and z we qualitatively 
expect a smoother field distribution for finite values of A compared to the 
hard-edge limit for A <^c 1. Note that the level of high spacial frequencies 
depends in a complicated way on Q, because the perfectly metallic waveg- 
uide effectively acts as a mirror. The value of A should be actually chosen to 
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Fig. 3. Comparison between hard-edge undulator case and high-frequency filtering. 
Function S(z') in the hard-edge case with A = (a) and in the case when filtering is 
applied with A = 0.1 (b). Their Fourier transforms !F{S} with respect toS, = C^'^ -I- C 
are plot in (c) for the hard edge case and in (d) when filtering is present. 

cut off high spatial frequencies that are outside the region of applicability 
of the resonance approximation, i.e. A ~ 1 /N^,. Thus, the correct value of A 
depends, case by case, on the type of setup considered. The Fourier trans- 
form of the function S(z') in Eq. (107) to be inserted into Eq. (108) and into 
Eq. (109) reads 



r [Sin + C)) = sine \l + C)\+^A exp f- f [Cf + cf 



X <cos 



i(er+c) 



sm 



^(cr+<:)]-fl[^(cr+<:) 
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. (110) 



where the imaginary error function erfi(-) is defined as 



IZ 

1 2 r 

erfi(z) - -erf(zz) = — — I exp [-f]dt . 
I iyn J 



(111) 



Fig. 3 presents a comparison between functions S and f{S} for A = 0.1 and 
A = 0, the latter case describing an undulator with hard edges. A filtering 
effect can be clearly seen, suppressing higher frequency components of !F{S}. 
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Hard-Edge Case 
X-polarization Y-polarization 




Fig. 4. Three-dimensional view of the virtual source in the free-space limit for A = 
(hard-edge case) and C = (perfect resonance) . 

Eq. (101) and Eq. (102) completely describe radiation from a planar undu- 
lator in the presence of a circular waveguide. The limit for Q » 1 can be 
taken as free-space limit, and allows one to compute the free-space case as 
a finite sum of modes, that is in an alternative way with respect to the cal- 
culation of integrals in [16]. Position z = corresponds to the virtual source 
discussed at the end of Section 2.4. The concept of laser-like radiation beam 
can in fact be naturally extended to the case a waveguide is present. Only, 
different boundary conditions have to be accounted for. Mathematically this 
means that propagation of the laser-like source must be performed with the 
proper Green's function, Eq. (74) multiplied by the numerical factor -licolc, 
rather than with the free-space Green's function, Eq. (83) (also multiplied 
by -2ico/c). As said before one must carefully select the number of modes 
used for computation, according to A: » Vo. In Fig. 4 we show a three- 
dimensional view of the virtual source in the free-space limit for A = and 
C = 0. We plot I = |Ep for both horizontal (x) and vertical (y) polarization 
components as a function of x and y, defined hy f = x Cx + y e^. 

We set C = (perfect resonance). In Fig. 5 we show a comparison between 
the analytic expression for the intensity distribution I at the virtual source 
position, that is obtained from a the normalized version of Eq. (30), and 
numerical expressions obtained through Eq. (101) for different values of 
A at Q = 1000 and Q = 100. A 2D plot is obtained by cutting the 3D 
intensity profile at y = (i.e. at (p = 0). Only the horizontal polarization 
component is important in this case. The importance of high-frequency 
filtering becomes clear in the case for Q = 100 and A = or A = 0.001. 
In general, the smaller the value of Q, the more filtering becomes relevant. 
Differences between results for A = and A = 0.001 should be taken as 
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Fig. 5. Intensity profiles of the virtual source (z = 0) at large values of Q = R^/iALip) 
for different values of A and C = (perfect resonance). This 2D plot is obtained 
cutting the 3D intensity profile at y = (i.e. at ^ = 0). The solid line is obtained 
with the help of Eq. (30). 

exemplification of the filtering process. However, they are irrelevant from 
a practical viewpoint. For example, if the number of undulator periods is 
of order 100, high frequency components for A = 0.001 or A = will not be 
distinguishable and, in fact, fall outside of the region of applicability of the 
resonance approximation. 

The behavior of the intensity profile is shown in Fig. 6 for A = 0.1 and in 
Fig. 7 for A = 0.001, where I = |Ep is plot for both the horizontal (Eq. (101)) 
and vertical (Eq. (102)) polarization components at z = 0. The number of 
modes used in all plots are in all cases much larger than Vo. Furthermore, 
we verified that results do not change by changing the number of modes 
used in the computation, provided that condition kmax » is fulfilled. 
Comparison of Fig. 6 and Fig. 7 clearly shows the effect of different edge 
dimensions. Note that in the free-space case, the horizontally polarized field 
is azimuthal symmetric. The situation becomes more complicated when the 
influence of tih.e waveguide begins to be important, i.e. for values of Q 
comparable with unity. In this case the field presents both horizontal and 
vertical polarization components, and azimuthal symmetry is lost, as is easy 
to see from figures and by inspection of Eq. (101) and Eq. (102). 

The knowledge of Eq. (101) and Eq. (102) solves the problem of charac- 



43 




Fig. 6. Intensity profiles of the virtual source (z = 0) at different values of 
Q = R^/(ALa,) for A = 0.1 and C = (perfect resonance). 

terizing planar undulator radiation in the presence of a circular vacuum 
pipe. Depending on the application one may use these equations in differ- 
ent ways. For example, one can propagate radiation in the presence of a 
vacuum pipe with a certain radius up to a given distance down the beam- 
line and subsequently propagate the electric field distribution in free-space 
with the help of Eq. (83) (or with the help of Fourier codes like ZEMAX 
|j22J). One may also account for changes of the waveguide radius by further 
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Fig. 7. Intensity profiles of the virtual source (z = 0) at different values of 
Q = K^/(ALj„) for A = 0.001 and C = (perfect resonance). 



solving the initial value problem for the field with the help of the proper 
Green's function, that will be proportional to Eq. (74). In principle one may 
even account for different geometries of the pipe, passing from a section 
of the beamline to another one, by investigating different eigenvalue prob- 
lems with respect to that for a circular pipe considered here. Many are the 
possibility of application of our theory. However, due to loss of azimuthal 
symmetry it is not easy to directly investigate the field distribution, because 
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plots of the field forcefully become three-dimensional, as we have just seen. 

It is thus worth to focus on finding figure of merits describing the influence 
of the pipe, that can give designers and scientists some measure of the 
influence of the waveguide in the stage of beamline design or planning of 
experiments. 

One figure of merit of interest is the ratio between the power density for 
a specific value of Q integrated over the waveguide cross-section and the 
angle-integrated power density in free space: 



-I 



£x(CJ) 




lim £x 



(112) 



When calculating the denominator in Eq. (112) we should introduce the 
same high-frequency cutoff. A, introduced for the numerator. This means 
that we substitute Eq. (108) and Eq. (109) (with S(z') defined in Eq. (107)) in 
Eq. (101) and Eq. (102), we take the limit for Q — > oo and we integrate in df± 
on a transverse plane at any fixed position z, the result being independent 
of z. In particular we choose to integrate the far-zone distribution and we 

make use of the definition 6 = fj_/z. Calling D the denominator in Eq. (112) 
we obtain: 
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The case for A = can be calculated analytically as: 
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(114) 



that also follows from a normalized version of Eq. (35). Numerical integra- 
tion gives D ^ 5.02 for A = 0.01 and D - 5.81 for A = 0.1, that are cases 
considered in this paper. It remains to calculate 
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(115) 
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In order to calculate the square modulus in Eq. (115) one has to compute the 
square modulus of an infinite sum in k. It can be shown that cross terms of 
this sum involving both TE and TM modes or different values of k vanish. 
One thus finds the result 



k=oo ( ^ 



Vq 



k=l 
k= 



vikr \ 



(116) 



Integrals in Eq. (116) can be calculated analytically. Further using the fact 
that Vik is the fc-th root of the Bessel function /i(-) we obtain 



^= ^ £ {i^^r + hiiiik) Mm) - uiiik)] + jiim)]] 

k=l 

+^ Z l'^^ [^'(^i*^^ + ^^M} ' (117) 

k=l 

where ^^iz) and ^l{z) are given in Eq. (108) and Eq. (109), while D is 
defined in Eq. (113). One may now study, for a given value A ~ lINyj, i.e. 
for a certain number of undulator periods, the dependence of W on the 
waveguide parameter Q or on the detuning parameter C. 

First we set C = and we plot, in Fig. 8, the function VV(Q) for different 
values of A ranging from to 0.1. It is seen that edge effects, that are very 
important in the intensity distribution, can be seen as a translation in Q 
when it comes to the figure of merit W(Q). Note that the dependence on 

Q in the coefficients \^''^ is of the form sin^ [iu2^/(4Q)] or sin^ [v^^/(4Q)] 

for C = 0, A = and z = 0. This expression is then multiplied by Q in Eq. 
(117). Thus W must go to zero as Q — > 0. However, for small values of 
Q, only a few waveguide modes are excited. One thus expects oscillations 
in W(Q). In particular when = n, i.e. at Q = 0.27 {^n = 1.84), we 

expect a value of W close to zero (not shown in Fig. 8). The first maximum 
of Qsin^ [jU^^/(4Q)], for values Q > 0.27 corresponds to Q = 0.73. The shift 
at higher values of Q in Fig. 8 is ascribed to contributions from other modes. 

Another possibility is to set a certain value for Q and to plot the function 
W(C) for different values of A. Since A ~ 1 /Na, it does not make sense to 
consider values C > 1 / A, because they are outside of the region of appli- 
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Fig. 8. Plot of VV as in Eq. (117) as a function of Q = B?-l{KLw) for C = (perfect 
resonance) at different values of A. 

cability of the resonance approximation. This explains the range of C in 
Fig. 9 and Fig. 10. In particular we fixed Q = 2.0 and we plot the case for 
A = 0.1 in Fig. 9 and values for A = 0.01 and A = in Fig. 10. It is interesting 
to remark the characteristic behavior in Fig. 10, where one can distinguish 
peaks at certain values of C. These peaks correspond to those values of C 
such that C^'^ + C = 0. At these values, the fc-th TE or TM mode is at reso- 
nance. Contrast of peaks becomes better as |C| increases, because different 
resonances are further away and peaks overlap less. For larger values of Q, 
as the number of excited modes increases, overlapping between different 
resonances hides the peaks, and one has the free-space limit (also shown for 
comparison in Fig. 10). Note that by definition of W, the integrated power 
density is normalized to the angle-integrated power density in free-space at 
resonance, i.e. at C = 0. This explains why VV = 1 at C = in Fig. 10. 

Finally, we may asses the influence of the waveguide on the output radiation 
by comparing the modulus of the vertical electric field for different values 
of Q with the modulus of the horizontally polarized field in free-space. In 
fact, as has been said before, the presence of boundary conditions different 
from those of free-space destroy the horizontal polarization as well as the 
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Fig. 9. Plot of W as in Eq. (117) as a function of C = 2nN^Aco/cor for Q = 2.0 at 
A = 0.1. 

azimuthal symmetry of the undulator radiation. Vertical and horizontal 
field are to be considered as functions of the distance from the z axis. Since 
the modulus of the vertical electric field is not azimuthal symmetric we take 
a cut in the direction where it is maximal, i.e. at (p = 7i/4, as it can be seen 
inspecting Eq. (102). As done before we limit our analysis for exemplification 
purposes to the case z = 0. We are thus interested in the functions Px and Py 
respectively defined as 



We studied Py as a function of f for different values of Q and different values 
of A. We set C — 0. Comparisons with (studied at Q = 1000) as a function 
of f at different values of A are presented in Fig. 11 for A = 0.1, Fig. 12 for 




(118) 



and 




(119) 
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Fig. 10. Plot of as in Eq. (117) as a function of C = InNia^cojcOr for Q = 2.0 at 
different values of A = 0.01 and A = 0. The far-field limit (Q — > oo) is also shown 
for comparison (dotted line). 

A = 0.01 and Fig. 13 for A = at different values of Q. Note that the maximal 
value of f is VO, and is due to the fact that the pipe poses a geometrical 
limit to the transverse region of interest. Figures underline the role of the 

vertical polarization component of the field for values of Q around unity. 
Thus, the level of the vertical polarization component of the field compared 
with that of the horizontal polarization gives a measure of the waveguide 
influence. 

As it is clear by inspection of Eq. (92) and Eq. (93), our method allows 
straightforwardly to study the evolution of the radiation pulse after the 
undulator, provided that the pipe section remains unvaried. In Fig. 14 we 
show the evolution of the intensity profile as a function of the normalized 
distance from the center of the undulator z. We take Q = 2.0, C = and 
A = 0.1. 

To conclude this Section we estimate the range of parameters of interest 
in the case of the infrared undulator beamline at FLASH. The wavelength 

range is between 50 f/m and 200 f/m. The radius of the vacuum pipe is 
R = 1.8 cm, the number of undulator periods N-a, = 9 with a period length 
Aj(, = 40 cm, yielding Ljy = 3.6 m. Note that the K parameter is always 
much larger than unity, increasing when the wavelength is increased by 
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Fig. 11. Comparison of Pj as defined in Eq. (118) (dotted line) and Py as defined 
in Eq. (119) (solid line) for A = 0.1 at different values of Q = R^/{KLu,) and C = 
(perfect resonance). Plots refer to the virtual source position (z = 0). 

ramping up the magnetic field in the undulator {K = 20 -v- 40 in the THz- 
gap operation range). This set of parameters means A ~ l/N-u, ~ 0.1 and 
3 ^ ^ ^ 10- As it can be seen in Fig. 6 and 11, for A = 0.1 in the short 
wavelength limit (Q ^ 10) the vacuum chamber has a small influence on 
the field distribution and polarization. A strong influence is visible, instead, 
in the long wavelength range (Q - 3). Finally, from Fig. 8 we see that the 
total power at perfect resonance (C = 0) is practically independent of the 
wavelength in our range of interest. 



7 Wall-resistance effects 



Up to now we studied the problem of undulator radiation in a waveguide 
assuming perfectly conducting walls, i.e. infinite conductivity o — > oo. 
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Fig. 12. Comparison of Px as defined in Eq. (118) (dotted line) and Py as defined 
in Eq. (119) (solid line) for A = 0.01 at different values of Q = R^/{XUv) and C = 
(perfect resonance). Plots refer to the virtual source position (z = 0). 

In this Section we follow the approach proposed in [|21|] to describe the 
case when the refractive index of the walls is still mainly defined by the 
conductivity o (i.e. we are still dealing with a metal), but o has a finite value. 
We will conclude that in practical situations of interest, the presence of wall- 
resistance effects introduce important changes to the theory for perfectly 
conductive walls. 

Let us start our investigations with Maxwell's equations in the time do- 
main, written in all generality for a medium having conductivity o and 
permittivity e: 



^ -i^ 4:no-t edE ^ -t IdH 

VxH= E + — r-, VxE = r- 

c c at c at 

V-H = 0, V-£ = 0. 



(120) 
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Fig. 13. Comparison of Px as defined in Eq. (118) (dotted line) and Py as defined in 
Eq. (119) (solid line) for A = (hard-edge) at different values of Q = R^KAL^p) and 
C = (perfect resonance). Plots refer to the virtual source position (z = 0). 

For a monochromatic wave of angular frequency co we have 

V X H = -i-{n'fE , VxE = i-H 
c c 

V-£ = 0, V-H = 0, (121) 

where 

I Ana 
n = \ e + 1 

V CO 

is the complex refractive index of the medium 
be written as 



(122) 

. A single equation for E can 
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Fig. 14. Intensity profiles at different values of z = zjL^ (z = at the undulator 
center) for Q = 2, C = (perfect resonance) and A = 0.1. 
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2 

V^E + {n'f^E = . (123) 



A solution of Eq. (123) is a plane wave specified by: 



E exp{-icot) + C.C. = Eg exp[z(fc' • r - cot)] + C.C. 

H^-{k'xE), (124) 

CO 

where k' = k'^Cx + k'yCy + k'^Cz is the (complex) wave vector inside the metal, 
with {k'f = {k'^f + {k'yf + {k'^f = {n')W/c^ a complex number. 

Let us consider a plane wave at the interface between vacuum (with unitary 

-> 

refraction index) and a medium described by Eq. (122), and let k be the wave 

-> 

vector of the incident wave in vacuum with k = \k\ = co/c. We indicate with 
subscripts t and n the tangential and normal components wave vectors to 
the metallic surface. Continuity of the tangential component of the wave 
vector on the boundary (i.e. k^ - kt) implies: 



4=1- 

As a result one obtains 




(126) 



In a metal, the refractive index is mainly defined by the conductivity o 
according to 

n' = Je + I ^ Ji . (127) 

\ CO y CO 

Since a » 1, also \n'\ » 1, and Eq. (126) yields 



K^^- (128) 

Thus, the propagation direction of the wave in the metal is almost per- 
pendicular to the surface, independently of the direction of the incident 
wave. This means that electric and magnetic field into the metal have only 
tangential components, Et and Ht. Moreover, these components must be 
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continuous on the surface, as follows from Maxwell's equations. Using Eq. 
(124) one obtains 



Eq. (129) can be considered as an approximate boundary condition. Note that 
this is valid not only for plane waves. In fact, any wave can be decomposed 
in terms of a linear superposition of plane waves. Since Eq. (129) is valid for 
each component, it must be valid for their linear superposition. Eq. (129) 
can also be extended for any shape of the boundary surface, provided that 
typical value of the curvature radius is much larger than the wavelength. 
Eq. (129) is named after Leontovich, who first derived it, and can be written 
in vector form as 



(nx£)| = -nx(nxH)| , 
Is Is 



(130) 



where the vector it, as defined as in Section 3 and shown in Fig. 2, is pointing 
inwards. Using Leontovich boundary condition gives a good approximation 
in the case of a metallic waveguide, and drastically simplifies the solution of 
the electrodynamical problem. Namely, the problem of mode excitation in a 
waveguide with resistive walls can be solved following the same approach 
in Section 3 and Section 4, where boundary conditions are now substituted 
by Eq. (130). 

Since we are applying the paraxial approximation, we assume that the field 
amplitude does not change much along the z direction within a wavelength 

and we are interested about the (main) transverse components of the field Ex 
and £y. Using Eq. (130) and Maxwell's equations, Eq. (121), it is possible to 
find boundary conditions involving the transverse field components only. 

-» 

First, using Eq. (121) in terms of slowly varying envelopes £, we can write 
Eq. (130) as 



(nx£j, =—^nx nx(vx£j 



(131) 



The transverse components of Eq. (131) only involve £z, while the compo- 

-» 

nent along Cz only involves £x. We may thus write 
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con' 



n X 



n X X Ex) I • 



(132) 



Eq. (132) will substitute the first boundary condition in (40). 

Second, from Eq. (130) we see that, on S, the projections of H and E on the 
plane orthogonal to it, i.e. Hf and £f, are orthogonal and their modulus are 
related by Eq. (129). As a result, we have £z = ^r^z = {l/n')Hf{ezXn).Now, 
Eq. (121) implies Vx -Ex = -(£l'/c)£z. Therefore Vx -Ex = -co/{n'c)Hf{ezXn) = 
-co/{n'c)H • (Cz X n). But Eq. (121) also yields H = -ic/coV x E. We then obtain 



i^^'H ^ii^""^)] -i^^x^-^n- (Vxl) xe 



(133) 



If we now calculate the expression in [•] brackets in Eq. (133) and we account 
for the fact that, in paraxial approximation, d^E^^y — {co/c)Ex^y » (?x,yEz we 
can re-write Eq. (133) in terms of slowly varying field amplitudes as 



(v.-l)|=i^(..l),, 



(134) 



that will substitute the second boundary condition in Eq. (40). 
Thus, the problem in (40) should be substituted with 



Ex(z,rx) 



= /(z,rx) 

X £x)| = -ic/icon') |n x n x ^Vx x £x j 
(Vx-£x)| ^ia)/{cn'){n-E\ . 



(135) 



We may now follow the same method used in Section 3 to solve (135). 
Namely, after consistent application of Laplace transformations, we find an 
eigenvalue problem analogous to (48): 



ViFy(r-l) + AjPjif'^) = 

(n X Fy)| = -iciicon') |n x [n x (Vx x F^)]| 

(Vx-F;)| =icol{cn')[fi-P)^ . 



(136) 
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As for the case of perfectly metallic walls, boundary conditions are homo- 
geneous, so that the domain of the Laplacian operator is the vector space of 
twice differentiable (square integrable) functions obeying boundary condi- 
tions in (136). However, the Laplacian operator defined in this way is not 
self-adjoint with respect to the inner product defined in Eq. (49). This is a 
result of the fact that vectors f y in Eq. (136), at the boundary, are not orthog- 
onal to the surface S. Then, eigenvalues are not real, nor eigenf unctions are 
orthogonal with respect to the inner product in Eq. (49). In general we do 
not know wether the spectrum is discrete, completeness is not granted and 
we cannot prove the existence of a set of eigenfunctions either. Yet, direct 



calculations in |21[] show that 

{F}Fi) = J Pj-Pid?^ = bji. (137) 



It follows that functions Fy form a bi-orthogonal set of eigenfunctions. Bi- 
orthogonality is often exploited in different problems (see e.g. [30, 31, 32]). 
In particular we will take advantage, without proving it, of completeness 
and discreteness of the bi-orthogonal set. This allows us to decompose G"^ 
as 

_ f.«f.„ 

G« = y ^. , (138) 

1^ Ljj liojpic -Aj ^ 



exactly as for the case of perfectly conducting walls. Note that now eigenval- 
ues and eigenfunctions are complex. In particular, we write the eigenvalues 
as 



Aj = A'. + iA'l . (139) 



We want to study the problem of wall-resistance in the framework of a per- 
turbation theory. If A° is an unperturbed eigenvalue (perfectly conductive 

walls) we require that \bAj\ = \Aj - A°\ « |A° - A°| for any value oiki^ j. 

In this case we may still formulate an eigenvalue problem for the two scalar 
(and complex) functions i/^J^ and i/'™ with complex eigenvalues Ay, using a 
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dejfinition analogous to Eq. (53). We get: 



ViV'J"'™(rx) + A]™ V^""'™(rx) = 
[it • V^xl'f + (e, X ii) • Vxi/^™]| = -ic (kff ^pfl{con') 

U Is 



(140) 



where 



i.TETM 



= 1 



(141) 



It should be noted that (140) is valid only within the framework of a pertur- 
bation theory. In fact, according to (140), functions xf^'™' separately obey 
Helmholtz's equation, but boundary conditions in (140) are not consistent 
with independency of ij^. ' required by Eq. (53). 

In the framework of a perturbation theory though, we may still use (140) 
and formally obtain the same tensor Green's function in Eq. (61). As said 
above, we want to work in the framework of perturbation theory. We base 
such theory on the small parameter wfi^n'k^c) «: 1. Within this theory, the 
problem in (140) can be written as 



Vii/.J^'™(rl) + \™ i/^f '™(rj = 

[n • V^i/;TE]| = |-fc (/cP)' ^]^f|{con') + ico/ [n'c (kffj [(el x H) ■ V^]' 
[{k™f ^™ I = ico/{cn') [ft • Vxi/'™]| , 



(142) 



where TE and TM modes are now decoupled. If we restrict our analysis to a 
circular waveguide radius R, we obtain the following expression for ipj ' 
at the first order in co/{{\n'kj_\c): 



TE 
mkl 



sin{m(p) ^ 
^ cos{m(p) 



(143) 



and 
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TM 
mkl ) 



{Vmk + ^Vmk) ^ 



sin(m(/)) 
^ cos{m(p) ^ 



(144) 



where -A^^'™ should now be calculated with the help of Eq. (141), while 
boundary conditions in (142) allow to obtain 



6|U 



mk 



ic ^l^ + co^m'RVicYj 



l^mk 



n'coR 



^Ik-^^ 



(145) 



and 



mk 



^mk 



icoR 
cn'v^ , 

mk 



(146) 



Note that Eq. (145) and Eq. (146) can be written in terms of our small 
parameter, since R(jol{c\n'\{v, = co/{c\n'k^l\), where the concept of trans- 
verse wave number is now applied to each mode k. Actually condition 
ijo/{{\n'k±\c) ^ 1 can be presented as \dvmk\ ^ 1 and \6[j.mk\ ^ 1- Thus, we 
may give an estimation of wall resistance effects retaining our main results, 
Eq. (101) and Eq. (102), and substituting [imk with [imk + ^l^mk and Vmk + bvmh 
In our case of interest, w = 1 so that 



fiik 



ic ^l^ + co^R^/{cYj 



n'coR 



(147) 



and 



^ = . (148) 

Addition of b\i\-k and 6vijc implies an exponential damping of each mode, 
due to the propagation factor exp[-zC]^''^z] in the expressions for J?l^'^(z), Eq. 
(103) and Eq. (104). We can estimate the damping exponent as 

^^.^ 1 L , 
^ \i\-^^\n'\R 



for TE modes, where we used {KIKf «; 1 and 

Dy = — - — z = — z (150) 

^ yf2\n'\R 
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for TM modes. Expressions for the damping exponents D^, and Dv, Eq. (149) 
and Eq. (150), are known in literature and maybe found, for example, in [33]. 
Note that the difference between D^, and Dy is only in the factor (/j^^ - 1)"^ 
in Eq. (149). D^, depends on the mode number, k, while is independent 
of it. 

Let us study a practical case for copper. We know that the skin depth for 
copper at wavelength A = 10 cm is 6s - 1.2 [im and 6s = VcA/ (In yfo), where 
a is the conductivity. It follows that for A = 200 /jm, that is in the range of 
interest for the FLASH infrared beamline we obtain 6s - 54 nm yielding 



V2A 

KI = f--840. (151) 

27i6s 

Let us consider R = 1.8 cm and Li„ = 3.6 m (corresponding to Q = 2.8), and 
substitute these values in Eq. (149) and Eq. (150). The perturbation parameter 
is smaller than unity: coR/{[iuc\n'\) - 0.37 (/Jn - 1.84), and coR/{viic\n'\) - 
0.18 (vii - 3.83). Therefore, we may apply the perturbation approach for 
estimations. It should be noted here that the accuracy of the first-order 
perturbation approach becomes better and better for higher modes, because 
[j-ik and vik increase while k increases. For the transverse magnetic modes 
we have a damping factor common to all modes Dy - O.V7z. For transverse 
electric modes we should specify the mode number. For the first mode 
{k = 1), and we obtain D^, - 0.07 z. For second mode (k = 2) we have 
/Ji2 - 5.3 and thus D^, - 0.006 z, that gives a negligible effect. All TM modes 
are visibly damped. At z = 1, i.e. half undulator length after the exit of the 
device, one has a damping factor exp(-0.17) - 0.84, i.e. modes are dumped 
by about 16%. TE modes instead, are left almost unaffected. 

We conclude that propagation through the pipe is strongly affected by wall- 
resistance effects even for relatively large values of Q. Note that this effect 
strongly depends on the material considered. For example steel has a skin 
depth three times thicker than copper. Strictly speaking we could not apply 
a perturbation approach, since now the perturbation parameter is larger 
than unity (although comparable with unity). However, if we still use the 
perturbation approach to get a rough estimation, we obtain that TM modes 
present a dumping exponent Dy that is about three times larger than copper. 
This means that TM modes are dumped by about 40%. These estimations 
allow us to formulate the following recommendation: the internal part of 
the vacuum pipe for the infrared undulator line at FLASH should be copper- 
coatec£il- 



^ We underline the fact that these remarks are only valid in the case of undulator 
radiation. In particular, as said before, we will treat edge-radiation setups in a 
separate work. 
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8 Conclusions 



In this article we presented the first exhaustive theory of undulator radia- 
tion within a waveguide and we exemplified it in the case of the infrared 
undulator beamline at FLASH. In the relatively simpler free-space case, 
paraxial Maxwell's equations can be solved in terms of a scalar paraxial 
Green's function. In that case, the solution for the field, written in reduced 
form within the resonance approximation is identical (aside for polariza- 
tion properties) for planar and helical undulator configurations. When a 
waveguide is present, the relation between current sources and electromag- 
netic field is more complicated. The analysis of the problem is performed 
by introducing a tensor Green's function technique, thus complicating the 
structure of equations that, in contrast to the free-space case, now depends 
on the undulator type and on the waveguide geometry. First we outlined a 
solution for a homogeneous waveguide with arbitrary cross-section. Then 
we specialized our consideration to the case of a circular waveguide, and 
we implemented the planar undulator cas^ within the applicability region 
of the resonance approximation. The electric field was found as a super- 
position of the waveguide modes, and was studied for different values of 
parameters. The main parameter involved in the problem is the waveg- 
uide parameter Q, that can be interpreted as the squared ratio between the 
waveguide radius and the radiation diffraction size, and is a purely geomet- 
rical parameter. When Q is comparable, or smaller than unity, waveguide 
effects become important, under the assumption of a perfect conductor. 
Moreover, we found that wall-resistance effects play a significant role even 
for relatively large values of Q. To minimize this effect, it is desirable to coat 
the internal part of the waveguide with copper. 
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We treated the helical case in Appendix A. 



62 



Appendix A: Helical undulator 

We define in all generality the horizontal and vertical velocity of an electron 
as a function of the longitudinal position z' along the undulator as 



Vx{z') - -cdssm{kwz') = — 2^ {exp[ik^uz'] - exp[-ik^uz']} 
Vy{z') = +c6sCos{k^z') = +^ {exp\iku,z'] + exp[-ik^z']} . 

(152) 

The + sign in the second equation in (152) indicates an electron rotating 
clockwise (- sign) or counterclockwise (+ sign) in the judgement of an 
observer located after the undulator and looking towards the device. The 
longitudinal Lorentz factor yz is now constant. We can write the phase factor 
in Eq. (62) as 



z 

I 



^^dz = -^-r = k^uZ + Cz . (153) 
2cyl 2Ayl 



The next step is to substitute Eq. (152) and Eq. (153) in Eq. (62) to calculate 
the field components. We can do that within the region of applicability of the 
resonance approximation, as in the case of a planar undulator. Within the 
region of applicability of the resonance approximation, condition (88), and 
neglecting negligible terms in Eq. (62) as specified in Section 5 we obtain the 
field components in cartesian coordinates for the case of a helical undulator: 



~ ^271^ J dz'(G? + /G«)^^,)^„exp[/Cz']. (154) 



By inspection of Eq. (74), under the resonance approximation only the terms 
with m = 1 survive. With the help of Eq. (74) and neglecting negligible terms 
we find for the horizontal field {a - 1) and for the vertical field {a - 2) 
respectively: 



^(z) Jo (vifc^) - Ji (vifc^) exp[+2i(p] 



+ 



(155) 
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and 



One can also write the expressions for the components in polar coordinates 
E,. and E^, that are related to and Ey through 

Er = E^ cos{(p) + Ey sm{(f)) ^^^^^ 
E^ = -Ej sm{(p) + Ey cos{(p) . 

We thus obtain 

„ oo 

E,.,ir,cp,z) = -i^Yj{^'^'^ 



exp[+i(p]\ (158) 



and 

„ oo 

£^,(r,c/),z) = T^2^tez) 
k=i ^ 

Note that both TM and TE modes are present in the expression for the field. 
It is interesting to compare this fact with what can be found in literature. 
As already remarked in Section 1, TM modes are neglected in [20]. For the 
readers' commodity we report here words from Section 9.3 of reference 
| [2o|l . There one may find: "It is easily seen that in the case of a spatially 
extended charge propagating in a helical wiggler, the TE modes (...) couple 
to the wiggler-induced motion, as the current components (...) By contrast, 
the TM modes (...) are driven by the uniform motion of the space-charge 
distribution in the cylindrical waveguide. Therefore, in the remainder of 
this derivation, we will focus on the TE modes". 

Neglecting TM modes is held by us as a misconception. If it was correct, 
also undulator radiation produced by planar undulators should exhibit 
only TE modes. In fact, as already discussed in Section 5, the field from a 
planar undulator can be seen as superposition of fields from two helical 



+ 



Jo (^I'c^) + h (vu^jexp[+2z(^] 



(156) 



Jo[h^lk^ + h(jilk^ 



l-lk 



exY>[+i(p\+ 



]o\Vlk^\ + ]2\Vik^ 



R 



RJ 



exp[+i(p] \ . (159) 
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trajectories. Without TM modes we could not even recover the free-space 
limit for R — > oo as studied in Section 5 for a planar undulator. In particular, 
we would loose both horizontal polarization and azimuthal symmetry of 
the field. 

Similarly, one can study the limit R — > oo for the helical case, following a 
procedure analogous to that in Section 5 for a planar undulator. In analogy 
with the planar case, one can show that if TM modes are neglected, both 
circular polarization and azimuthal symmetry of the field are lost. 

Neglecting TM modes in f2d\ is the result of a more general misconception 
introduced when Maxwell's equations are solved with the Green's func- 
tion approach. Namely, the misconception is introduced in the eigenmode 
decomposition of the four-current (see paragraph 5.6.4 of 120!]), where TM 
modes are not coupled with the transverse current density as they must be. 
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